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Current.  Status _of_the  Developroent  o£  an  Ignition 
Transient  Model  for  Solid  Rocket  Motors 


Gary  D.  Luke 
Aerojet  ASRM  Division 
Sacramento,  California 
and 

Harry  A.  Dwyer 
University  of  California 
Davis,  California 


An  ignition  transient  model  for  solid  rocket  motors  is  currently  being 
developed  jointly  by  the  Aerojet  ASRM  Division  and  the  University  of  Califor- 
nia at  Davis.  Though  the  CFD  code  will  be  general  enough  to  predict  the  start 
transient  for  most  solid  rocket  motors,  the  particular  motive  for  the  develop- 
ment of  this  code  stems  from  the  desire  to  analyze  the  flow  field  within  the 
Advanced  Solid  Rocket  Motor  (ASRM)  for  the  Space  Shuttle  with  all  its  geomet- 
ric complexity.  The  star  grain  configuration  in  the  head  end  of  the  ASRM 
coupled  with  the  multiport  igniter  creates  a formidable  problem  which  can  only 
be  modeled  accurately  using  a three  dimensional  Navier-Stokes  code  if  one 
wishes  to  preserve  both  volume  and  burning  surface  area  as  a function  of  axial 
position  down  the  bore.  The  actual  physical  geometry  is  crucial  in  modeling 
the  multiple  wave  interactions  occurring  within  the  combustion  chamber  as  well 
as  in  predicting  the  correct  amount  of  mass,  momentum,  and  energy  injected  as 
a function  of  time  and  space  at  the  propellant  surface. 

The  primary  objectives  of  the  CFD  code  are  to  calculate  the  pressure  rise 
rate,  the  thrust  rise  rate,  and  the  maximum  chamber  pressure  which  occurs 
during  the  first  second  of  the  ASRM  action  time.  And,  more  specifically,  to 
determine  the  relative  difference  between  the  ignition  transients  produced 
using  a single  port  igniter  and  that  produced  by  a multiport  igniter. 

An  implicit,  three  dimensional,  time  accurate,  finite  volume  method  is 
being  developed  to  solve  this  problem.  Current  plans  are  to  use  an  ADI 
technique,  with  replacement,  to  solve  the  set  of  analytically  linearized  full 
Navier-Stokes  equations.  However,  to  gain  an  understanding  of  which  features 
of  the  code  need  more  attention  than  others,  and  to  provide  an  inexpensive, 
quick  tool  for  studying  the  time  accuracy  of  the  selected  solution  algorithm, 
a one  dimensional  version  of  the  code  was  written  first.  The  various  features 
of  the  one  dimensional  code  were  validated  by  comparing  the  numerical  results 
to  analytical  results  for  problems  where  exact  solutions  were  available.  Then 
the  one  dimensional  code  was  used  to  perform  some  sensitivity  studies  to  help 
develop  an  understanding  of  the  complex  wave  interactions  occurring  within  the 
solid  rocket  motor. 

This  presentation  will  discuss  the  results  obtained  from  the  one  dimen- 
sional code  and  the  plans  for  the  development  of  the  3-D  code. 
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CURRENT  STATUS  OF  THE  DEVELOPMENT  OF  AN 
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Single  (RSRM)  and  Multiport  (ASRM) 
Igniters 


DEFINITION  OF  PROBLEM 

NEED  TO  ACCURATELY  PREDICT: 
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Details  of  Flow  Field  in  Head  End 
Where  Fins  and  Igniter  are  Located 
for  RSRM  and  ASRM 


THREE  PART  MODELING  APPROACH 
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PURPOSE  OF  1-D  CODE 
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LOCXHEED  • AEROJET  • ROST 


DESCRIPTION  OF  NUMERICAL  MODEL 
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System  of  Equations  Solved  Using  Alternating 
Direction  Implicit  (ADI)  with  Replacement 
(Essentially  a Direct  Solver  for  1 -D  Case) 


Governing  Equations  in  Integral  Form 
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LOCKHEED  • AEROJET  • RUST 


Linearized  Form  of  Discrete  Equations 


Analytically  Linearized  Using  Newton’s  Method 


Artificial  Dissipation 

2nd  Order  Terms  Required  to  Dampen  Out 
Oscillations  Occurring  at  Steep  Gradients 
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FEATURES  OF  1-D  CODE 
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Mass  Injection  Through 
Side  Walls 
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- Impermeable 

- Specified  Mass/Energy  Injection  Rate 

- Burning  Propellant  Surface 


VALIDATION  BY  COMPARISON 
TO  EXACT  1-D  SOLUTIONS 
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Shock  Tube  Problem 


Flow  Through  Constant  Area  Duct  With  Friction 
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LOCKHEED  • AEROJET  • RUST 


I-LUW  THROUGH  CONSTANT  AREA  DUCT  WITH  FRICTION 
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Numerical  results  from  AJITC1D  Code  using  200  cells 
Analytical  results  interpolated  from  Fan  no  Flow  Tables 


f-LUW  IHHOUGH  CONSTANT  AREA  DUCT  WITH  FRICTION 
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Numerical  results  from  AJITC1D  Code  using  200  cells 
Analytical  results  interpolated  from  Fanno  Flow  Tables 


ISENTROPIC  FLOW  THROUGH  A NOZZLE 

NOZZLE  CONFIGURATION 
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Numerical  results  from  AJITC1D  Code  using  145  ceils 
Analytical  results  from  l-D  Isentropic  Tables 


ISENTROPIC  FLOW  THROUGH  A NOZZLE 

PRESSURE  vs  AXIAL  POSITION 


•Bred  ‘aanssaud 
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Numerical  results  from  AJITC1D  Code  using  145  cells 
Analytical  results  from  1-D  Isentroplc  Tables 


ISENTROPIC  FLOW  THROUGH  A NOZZLE 

TEMPERATURE  vs  AXIAL  POSITION 


U68p  •3bniVU3dW31 


743 


Numerical  results  from  AJITC1D  Code  using  145  cells 
Analytical  results  from  1-D  Isentropic  Tables 


Radial  Mass  Injection  Into  Constant  Area  Duct 
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RADIAL  MASS  INJECTION  INTO  CONSTANT  AREA  DUCT 
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Numerical  results  from  AJITC1D  Code  using  200  cells 
Analytical  results  from  i-D  formulas  in  Shapiro's 


nnuiML  MASS  INJECTION  INTO  CONSTANT  AREA  DUCT 
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Numerical  resultsfrom  AJITC1D  Codeusing  200  cells 
Analytical  results  from  1-D  formulas  in  Shapiro's 


ONE  DIMENSIONAL  SHOCK  TUBE  PROBLEM 

PRESSURE  vs  AXIAL  POSITION  PRIOR  TO 
SHOCK  REFLECTION,  t = 0.5  milliseconds 
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DIAPHRAGM  ORIGINALLY  LOCATED  AT  X/L  = 0.75 
LEFT:  P = 400  psia,  T = 530  degR,  V = 0 ft/sec 
RIGHT:  P = 100  psia,  T = 530  degR,  V = 0 ft/sec 


ONE  DIMENSIONAL  SHOCK  TUBE  PROBLEM 

PRESSURE  vs  AXIAL  POSITION  AS  SHOCK 
CONTACTS  WALL,  t = 1.0  milliseconds 
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DIAPHRAGM  ORIGINALLY  LOCATED  AT  X/L  = 0.75  f^VblCf  Y' 

LEFT:  P = 400  psia,  T = 530  degR,  V = 0 ft/sec 
RIGHT:  P = 100  psia,  T = 530  degR,  V = 0 ft/sec 


ONE  DIMENSIONAL  SHOCK  TUBE  PROBLEM 

PRESSURE  vs  AXIAL  POSITION  NEAR 
STEADY  STATE,  t = 835  milliseconds 
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AXIAL  POSITION,  (in.) 

DIAPHRAGM  ORIGINALLY  LOCATED  AT  X/L  = 0.75 
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ONE  DIMENSIONAL  SHOCK  TUBE  PROBLEM 

TEMPERATURE  vs  AXIAL  POSITION  PRIOR  TO 
SHOCK  REFLECTION,  t = 0.5  milliseconds 


RIGHT:  P = 100  psia,  T = 530  degR,  V = 0 It/sec 
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ONE  DIMENSIONAL  SHOCK  TUBE  PROBLEM 

TEMPERATURE  vs  AXIAL  POSITION  AS  SHOCK 
CONTACTS  WALL,  t = 1.0  milliseconds 


AXIAL  POSITION,  (in.) 

DIAPHRAGM  ORIGINALLY  LOCATED  AT  X/L  = 0.75 
LEFT:  P = 400  psia,  T = 530  degR,  V = 0 ft/sec 
RIGHT:  P = 100  psia,  T = 530  degR,  V = 0 ft/sec 
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ONE  DIMENSIONAL  SHOCK  TUBE  PROBLEM 

TEMPERATURE  vs  AXIAL  POSITION  NEAR 
STEADY  STATE,  t = 835  milliseconds 


PRELIMINARY  SENSITIVITY  STUDIES 

MAXIMUM  HEAD  END  dP/dt  AS  A FUNCTION  OF 
DUCT  LENGTH  AND  NOZZLE  CONTRACTION  RATIO 
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PRELIMINARY  RESULTS  FROM  THE  AJITC1 D CODE 

VERIFICATION  OF  SHOCK  STEEPENING  AS  A 
CONSEQUENCE  OF  AN  INCREASING  MASS  FLOW  RATE 
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IGNITER  PRESSURE  RISE  RATE  = 75,000  psi/sec 


PRELIMINARY  RESULTS  FROM  THE  AJITC1 D CODE 

VERIFICATION  OF  SHOCK  STEEPENING  AS  A 
CONSEQUENCE  OF  AN  INCREASING  MASS  FLOW  RATE 
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PRELIMINARY  RESULTS  FROM  THE  AJITC1 D CODE 

MAXIMUM  dP/dt  AS  A FUNCTION  OF  THE  TOTAL 
DUCT  LENGTH  AND  THE  IGNITER  MASS  FLOW  RATE 
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PRELIMINARY  RESULTS  FROM  THE  AJITC1 D CODE 

MAXIMUM  dP/dt  AS  A FUNCTION  OF  THE  TOTAL 
DUCT  LENGTH  AND  THE  IGNITER  MASS  FLOW  RATE 
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PRELIMINARY  RESULTS  FROM  THE  AJITC1 D CODE 

MAXIMUM  dP/dt  AS  A FUNCTION  OF  THE 
DOWNSTREAM  BOUNDARY  CONDITION  SPECIFIED 
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RESULTS  OBTAINED  USING  AN  INCREASING 
IGNITER  FLOW  RATE  (75,000  PSI/SEC) 


CURRENT  STATUS  OF  1-D  CODE 
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ASRM  DIVISION 


FUTURE  PLANS  FOR  1-D  CODE 
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ASRM  DIVISION 


PLANS  FOR  3-D  IGNITION  TRANSIENT  MODEL 
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SRMAFTE  Facility  Checkout  Model  Flow  Field  Analysis 

Richard  A.  Dill,  ERC  Incorporated 
Harold  R.  Whitesides,  ERC  Incorporated 

Abstract 

The  Solid  Rocket  Motor  Air  Flow  Equipment  (SRMAFTE)  facility  was  constructed  for  the  purpose  of 
evaluating  the  internal  propellant,  insulation,  and  nozzle  configurations  of  solid  propellant  rocket  motor  designs. 
This  makes  the  characterization  of  the  facility  internal  flow  field  very  important  in  assuring  that  no  facility 
induced  flow  field  features  exist  which  would  corrupt  the  model  related  measurements.  In  order  to  verify  the 
design  and  operation  of  the  facility,  a three-dimensional  computational  flow  field  analysis  was  performed  on  the 
facility  checkout  model  setup. 

The  facility  checkout  model  entails  a straight  constant  diameter  pipe  in  place  of  a specific  solid  propellant 
rocket  motor  internal  component.  This  configuration  was  to  provide  a simple  internal  flow  field  for  evaluation  of 
any  facility  induced  effects. 

One-dimensional  estimates  of  the  checkout  model  flow  field  were  available  for  comparison  to  the 
measurement  data  collected  for  the  checkout  model  but  the  CFD  results  provided  a comparative  estimate  in 
regions  where  one-dimensional  estimates  were  not  valid. 

Since  the  facility  was  too  large  and  complex  to  perform  a complete  three-dimensional  analysis  from  end 
to  end,  the  facility  was  divided  into  three  major  zones  for  analysis.  1)  The  header  pipes,  metering  nozzle,  nozzle 
exit,  and  diffuser.  2)  The  adapter  chamber,  transition,  checkout  model  section,  and  checkout  model  nozzle.  3) 
The  model  nozzle  exit,  diffuser  and  muffler.  The  three-dimensional  numerical  calculation  of  the  flow  field  was 
performed  by  Fluent/BFC.  This  code  solves  the  full  Navier-Stokes  equations  of  fluid  flow  cast  in  a staggered 
grid  formulation.  The  SIMPLER  numerical  algorithm  is  used  in  the  solution  process.  The  code  utilizes  wall 
functions  instead  of  physically  resolving  the  boundary  layer  and  the  standard  k-e  turbulence  model  is  used  to  close 
the  system  fluid  dynamic  equations. 

The  checkout  model  measurement  data,  one-dimensional  and  three-dimensional  estimates  were  compared 
and  the  design  and  proper  operation  of  the  facility  was  verified.  The  proper  operation  of  the  metering  nozzles, 
adapter  chamber  transition,  model  nozzle  and  diffuser  were  verified.  The  one-dimensional  and  three-dimensional 
flow  field  estimates  along  with  the  available  measurement  data  are  compared  in  this  presentation. 
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STATIC  PRESSURE  (psia)  600 


MPUTATIONAL  GRID  FOR  MANIFOLD  SYSTEM  AND 
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OUNDARY  CONDITIONS  USED  FOR  THE 

DEL  ADAPTER.  SPOOL  PIECES  AND  MODEL  NOZZLE 
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Axial  Direction  Grid 
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-D  FHFCKOUT  MODEL  FLOW  FIELD  (21X30X170)  GRID 


Pressure  Ratio  and  Mach  Number  Profiles  at  Model  Nozzle 

Throat  Plane 


Normalized  Radial  Distance,  r/Rwall 


UNDARY  CONDITIONS  USED  FOR  THE  SRMAFTE  DIFFUSER 
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Model  Aft  Total  Pressure/ Ambient  Pressure 
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A Comparative  Study  of  the  Effects  of  Inhibitor  Stub  Length  on 
Solid  Rocket  Motor  Combustion  Chamber  Pressure  Oscillations: 

RSRM  at  T=80  Seconds,  Preliminary  Results 

D.  Chasman,  D.  Burnette,  J.  Holt 
Rockwell  International,  Space  Systems  Division 
Huntsville,  Alabama  35806 
R.  Farr 

NASA  Science  and  Engineering  Laboratory 
George  C.  Marshall  Space  Flight  Center 

ABSTRACT 

Results  from  a continuing,  time-accurate  computational  study  of  the 
combustion  gas  flow  inside  the  Space  Shuttle  Redesigned  Solid  Rocket  Motor 
(RSRM)  are  presented.  These  CFD  analyses  duplicate  unsteady  flow  effects  which 
interact  in  the  RSRM  to  produce  pressure  oscillations,  and  resulting  thrust 
oscillations,  at  nominally  15,  30  and  45  Hz.  Results  of  Navier-Stokes 
computations  made  at  mean  pressure  and  flow  conditions  corresponding  to  80 
seconds  after  motor  ignition  both  with  and  without  a protruding,  rigid  inhibitor  at 
the  forward  joint  cavity  are  presented  here. 

Previous  studies  by  the  authors  have  demonstrated  that  combustion  chamber  pressure 
oscillations  in  the  RSRM  are  generated  by  flow/acoustic  interactions  which  occur  at  the  three  field 
joint  cavities  [1,2,3].  Edge-tone,  Hole-tone  and  Organ  pipe-tone  are  all  acoustic  sources  that  play 
part  in  these  interactions[4,5,6].  By  constructive  interference  proccesses,  the  different  acoustic 
sources  interact,  amplifying  the  pressure  oscillation  level  at  some  instant  during  the  RSRM  bum 
(i.e.  T+80  sec).  This  behavior  is  representitive  of  an  aerodynamic  whistle  of  Class  III  [7,8], 

However,  the  question  remained  as  to  whether  the  cavities  alone  are  the  main  acoustic 
generators  or  whether  protruding  inhibitors  dominate  the  system.  With  this  in  mind,  two 
simulations  have  been  conducted.  The  first  simulation  represent  a full  scale  model  of  the  RSRM 
with  a rigid  inhibitor  perturbing  the  flow  at  the  forward  joint,  while  the  second  simulation  was 
conducted  without  an  inhibitor.  All  other  flow  conditions  and  grain  geometries  were  kept  constant 
for  both  simulations. 

Fig.  1 shows  a comparison  of  results  between  the  two  simulations.  A full-length  density 
contour  plot  of  the  simulation  with  forward  inhibitor  is  shown  in  Fig.  la.,  while  that  without  the 
inhibitor  is  shown  in  Fig.  Id..  In  both  plots  the  shear  layer  which  developed  from  the  burning 
surfaces  is  apparent.  High  density  values,  corresponding  to  areas  of  high  pressure  amplitude, 
appear  in  red. 

Details  of  the  forward  field  joint  area  are  shown  for  both  simulations  in  Figs.  lb.  and  le.. 
Both  show  intense,  but  different,  shear  flow  activity  indicative  of  vortex  dynamics.  Streamlines 
shown  in  Fig.  lc.  and  Fig.  If.  single  out  individual  vortices  and  further  illustrate  the  difference 
between  the  two  simulations. 

Fig.  2 illustrates  a comparison  of  pressure  data  time  histories  for  three  points  common  to 
both  simulations.  It  can  be  seen  that  30  Hz  oscillations  are  evident  in  the  midsection  of  the  chamber, 
while  15  Hz  is  found  at  both  the  head  and  aft  ends,  with  the  head  end  being  180  degrees  out  of 
phase  with  the  aft  end.  These  results  indicate  classic  organ  pipe  acoustics  are  found  in  both 
simulations.  However,  there  is  a marked  difference  in  the  peak-to-peak  amplitude  of  these  pressure 
oscillations.  These  preliminary  results  show  pressure  amplitude  in  the  case  without  the  inhibitor 
are  about  14%  of  total  chamber  pressure,  while  those  in  the  case  with  the  inhibitor  are  only  10%  of 
the  total  pressure. 

While  calculated  peak-to-peak  pressure  amplitude  values  are  significally  higher  than  levels 
measured  during  actual  firing  and  flight,  we  feel  these  results  can  nevertheless  be  used  for 
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comparative  analyses  of  the  effects  of  inhibitor  stub  height  on  RSRM  combustion  chamber  acoustic 
pressure  amplitudes.  Specifically,  these  findings  demonstrate  that  the  inhibitors  alone  are  not  the 
dominant  factor  in  amplifying  combustion  chamber  pressure  oscillations,  but  rather  are  included  in 
secondary  acoustic/flow  interactions  occuring  at  the  three  field  joint  cavities.  Our  results  indicate 
that  inhibitors,  when  present,  actually  act  to  damp  such  oscillations. 
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RSRM  PC  OSCILLATION 


b)  Density  Detail  c)  Streamline  Detail 

With  Forward  Inhibitor  With  Forward  Inhibitor 


e)  Density  Detail 
Without  Forward  Inhibitor 


f)  Streamline  Detail 
Without  Forward  Inhibitor 


Figure  1.  Simulations  With  and  Without  Forward  Inhibitor 
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Overview  of  the  relevant  CFD  work 
at  Thiokol  Corporation 

Pawel  Chwalowski  & Hai-Tien  Loh 
Thiokol  Corporation,  M/S  L63,  P.0.  Box  707, 

Brigham  City,  Utah  84302-0707 


The  use  of  computational  fluid  dynamics  (CFD)  in  supporting 
the  rocket  propulsion  designs  at  Thiokol  Corporation  has 
continuously  increased  in  the  past  few  years.  An  in-house 
developed  proprietary  advanced  CFD  code  called  SHARP®  is  a 
primary  tool  for  many  flow  simulations  and  design  analyses. 
The  SHARP  code  is  a time  dependent,  two-dimensional  (2-D) 
axisymmetric  numerical  solution  technique  for  the  compressible 
Navier-Stokes  equations.  The  solution  technique  in  SHARP  uses 
a vectorizable  implicit,  second  order  accurate  in  time  and 
space,  finite-volume  scheme  based  on  the  following:  1)  an 
upwind  flux-difference  splitting  of  a Roe-type  approximated 
Riemann  solver,  2)  Van  Leer's  flux  vector  splitting,  and  3)  a 
fourth  order  artificial  dissipation  scheme  with  a 
preconditioning  to  accelerate  the  flow  solution.  Turbulence 
is  simulated  by  an  algebraic  model,  and  ultimately  the  k-e 
model.  Some  other  capabilities  of  the  code  are  2-D  two-phase 
Lagrangian  particle  tracking  and  cell  blockages.  Extensive 
development  and  testing  has  been  conducted  on  the 
three-dimensional  (3-D)  version  of  the  code  with  flow, 
combustion,  and  turbulence  interactions. 

The  SHARP  code  has  been  applied  in  many  areas  of  the  solid 
rocket  motor  ( SRM ) design  involving  internal  and  external  flow 
analysis.  However,  the  internal  flow  analysis  inside  the 
motor  and  in  the  nozzle  region  are  the  most  frequent. 
Usually,  the  results  from  these  CFD  analyses  become  the 
boundary  conditions  for  thermal  and  structural  computations. 
In  the  case  of  the  internal  nozzle  flow  calculations,  SHARP 
computes  the  convective  heat  transfer  coefficients  and 
temperature  distribution  along  the  nozzle  wall  for  the  thermal 
erosion  predictions,  and  the  pressure  distribution  for  the 
structural  predictions.  The  pressure  loads  on  the  propellant 
grain  surfaces  inside  the  SRM  obtained  from  SHARP  are  used  as 
a boundary  conditions  to  predict  propellant  grain  deformation 
and  displacement. 

The  2-D  two-phase  Lagrangian  particle  tracking  gives  the 
ability  to  predict  solid  particle  impingement  on  the  exit 
cone.  Also,  SHARP  prediction  of  the  slag  accumulation  in  the 
aft  dome  region  of  the  SRM  agrees  with  the  actual  static  test 
data . 

The  emphasis  in  the  presentation  will  be  put  on  the  specific 
applications  of  SHARP  in  SRM  design. 


® SHock  wave  And  Recirculation  Program  is  a copyrighted 
acronym  owned  by  Thiokol  Corporation. 
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3D  version  in  development  and  testing 


SHARP  2D  DESCRIPTION 
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APPLICATIONS  IN  SRM  DESIGN 
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SHARP  Results 
Internal  Flow  Modeling 
Nozzle  Flow 

Pressure  Contour  (N/m2) 
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SHARP  Results 
Internal  Flow  Modeling 
Nozzle  Flow 
Mach  Contour 
Temperature  Contour  (K) 


INFORMATION  ON  THIS  PAGE  WAS  PREPARED  TO  SUPPORT  AN  ORAL  PRESENTATION 
AND  CANNOT  BE  CONSIDERED  COMPLETE  WITHOUT  Tl€  ORAL  DISCUSSION 


APPLICATIONS  IN  SRM  DESIGN  (cont. 
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Internal  Flow  Modeling 
Shuttle  SRM  Internal  Flow 
Velocity  Vector  Contour  (ft/sec) 
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Submitted  for  the  CFD  Workshop  - 1992 

A Status  of  the  Activities  of  the 
NASA/Marshall  Space  Flight  Center 
Combustion  Devices  Technology  Team 

Kevin  Tucker 


The  Consortium  for  Computational  Fluid  Dynamics  (CFD) 
Applications  in  Propulsion  Technology  was  established  to  focus 
CFD  applications  in  propulsion.  Specific  areas  of  effort 
include  developing  the  CFD  technology  required  to  address  rocket 
propulsion  issues,  validating  the  technology,  and  applying  the 
validated  technology  to  design  problems;  all  under  peer  review 
by  experts  in  the  field. 

The  Combustion  Devices  Technology  Team  was  formed  to  implement 
the  above  objectives  in  the  broad  area  of  combustion-driven 
flows.  In  an  effort  to  bring  CFD  to  bear  in  the  design 
environment,  the  team  has  focused  its  efforts  on  the  Space 
Transportation  Main  Engine  nozzle.  The  main  emphasis  has  been 
on  the  film  cooling  scheme  used  to  cool  the  nozzle  wall. 
Benchmark  problems  have  been  chosen  to  validate  CFD  film  cooling 
capabilities.  CFD  simulations  of  the  subscale  nozzle  (to  be 
tested  8/92)  have  been  made.  Also,  CFD  predictions  of  the  base 
flow  resulting  from  this  type  of  nozzle  have  been  made.  A 
status  of  these  calculations  will  be  presented  along  with  future 
plans. 


National  Aeronautics  and 
Space  Administration 


KEVIN  TUCKER 
APRIL  29, 1992 
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Remtech 

University  of  Tennessee  Space  institute  (UTSI) 
Pennsylvania  State  University 
The  University  of  Alabama  (UA) 

The  University  of  Alabama  in  Huntsville  (UAH) 
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- High  aspect  ratio  coolant  channels 

- Compared  general  technology  issues  with  STME  TCA  concerns 

- Chose  STME  film/dump  cooled  nozzle  for  team  focus 

- Identified  validation/verification  requirements  for  supersonic  film  cooling 
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- Validate  CFD  codes  for  supersonic  film  cooling 

• Use  validated  codes  as  design  tool  for  subscale  (40K)  nozzle 

- Verify  film  cooling  performance  in  subscale  nozzle 

- Use  validated  codes  in  design  of  full  scale  nozzle 
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- Demonstration  of  CFD  as  a nozzle  design  tool  successful 

- Subscale  secondary  injector  redesigned  using  CFD 

- Confidence  gained  in  film  cooling  early  in  subscale  design 

- Large  film  cooling  analytical/test  database  being  developed 
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CORE/SECONDARY  COOLANT  FLOW  INTERACTION 

Interaction  Pressure  Contours  Indicate  a Step— Induced  Shock 
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CORE/SECONDARY  COOLANT  FLOW  INTERACTION 

Velocity  Vectors  / H 2 Mass  Fraction  Contours  Show  Vortex  Mixini> 


CORE/SECONDARY  COOLANT  FLOW  INTERACTION 

Interaction  Temperature  Contours  Indicate  Entrainment  of  Core  Gases 
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CORE/SECONDARY  COOLANT  FLOW  INTERACTION 

II 2 Muss  Fraction  Contours  Indicate  Vortex  and  Diffusion  Mixing 
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NASA 

National  Aeronautics  and  Combustion-Driven  Flow  Analysis  Technology 

Space  Administration 
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Impact 

- Results  of  "demonstration"  calculations  delivered  to  ED64  for  heat  flux  calculations 

- The  3D  analyses  are  pushing  computer  limits;  job  turnaround  is  slow 
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VIEW  LOOKING  FWD 
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space  SI"  Combustion-Driven  Flow  Analysis  Technology 
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0.0  0.4  0.0  1.2  I.B  2.0  2.1  2.B  3.2  3.6 

distance  from  center  of  heat  shield.  Inches 
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CFD  contributing  to  STME  design  improvements 


N92-32251 


CFD  Analysis  of  the  STME  Nozzle  Flowfield 
Anantha  Krishnan 

CFD  Research  Corporation,  Huntsville,  AL 

and 

Kevin  Tucker 

NASA  MSFC,  Huntsville,  AL 


The  Space  Transporter  Main  Engine  (STME)  uses  a gas  generator  cycle  to  cool  the 
nozzle  wall  by  a film-dump  of  the  turbine  exhaust.  The  ability  to  cool  the  skirt  is  a 
key  concern  in  the  design  of  the  STME.  CFD  calculations  were  undertaken  to  predict 
the  film  cooling  effectiveness  and  performance  sensitivities  for  various  design 
configurations,  operating  points  and  inlet  conditions.  The  results  in  this  study  were 
obtained  for  the  subscale  nozzle.  The  computations  were  performed  using 
REFLEQS. 

The  computational  analysis  showed  that  a chemical  equilibrium  model  was 
necessary  to  obtain  correct  predictions  of  the  specific  impulse.  The  frozen  composi- 
tion model  underpredicts  the  ISP  by  about  6%.  It  was  also  observed  that  the  coolant 
film  was  successful  in  maintaining  the  nozzle  wall  well  below  the  stagnation 
temperature  of  the  core  flow.  The  effect  of  the  coolant  flow  on  the  performance  of 
the  engine  was  found  to  be  negligible.  The  computed  heat  fluxes  at  the  wall  were  in 
good  agreement  with  the  empirical  data  obtained  by  Pratt  & Whitney. 

Further  test  data  from  Pratt  & Whitney  are  forthcoming  for  the  subscale  nozzle. 
Calculations  will  be  performed  to  determine  cooling  efficiencies  and  nozzle 
performance  over  a range  of  conditions  and  the  model  predictions  will  be  compared 
with  experimental  data. 
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Tenth  Annual  CFD  Workshop 
NASA  - MSFC 
April  28-30, 1992 
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The  Calculations  were  Done  for  the  Subscale  Nozzle. 
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CFD  METHODOLOGY 
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Validated  for  Large  Number  of  Benchmark  Problems 
(REFLEQS  Validation  Manual) 
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CASES  CONSIDERED 
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COMPUTATIONAL  RESULTS  FOR  THE 
SUB-SCALE  NOZZLE  FILM  COOLING 
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- P&W  Analysis  for  Primary  Injector 

- P&W  Analysis  for  Secondary  Injector 
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The  Nozzle  Performance  is  Relatively  Insensitive 
to  the  Coolant  Injection 


N92-32252 

NLS  Nozzle  Base  Flow  Characteristics 

J.  J.  Erhart 
Pratt  & Whitney 
West  Palm  Beach,  FL 

ABSTRACT 

The  flow  characteristics  of  the  NLS  nozzle  base  area  need  to  be 
determined  in  order  for  heat  transfer  rates  to  be  estimated.  The 
objective  of  this  work  is  to  calculate  these  flow  characteristics 
using  CFD . A Full  Navier-Stokes  code  in  an  axisymmetric  mode  using  a 
k-e  turbulent  model  with  wall  functions  is  applied.  Calculations  were 
completed  at  an  altitude  of  3,250  and  80,000  feet  in  the  flight 
trajectory.  The  results  show  flow  features  which  can  affect  vehicle 
design.  Calibration  of  a 3-D  case  with  data  is  underway. 
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MOTIVATION 

Understand  Base  Flow  Phenomena 
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GEOMETRY  & FLOW  BOUNDARY  CONDITIONS 
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RESULTS 

H2  - Low  Altitude 
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RESULTS 

Mach  Number—  Low  Altitude 
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SUMMARY 


N92- 32253 

HEAT  TRANSFER  IN  ROCKET  ENGINE 
COMBUSTION  CHAMBERS  AND  NOZZLES 

P.  G.  Anderson*,  Y.S.  Chenf,  and  R.C.  Farmer* 

Abstract 


Complexities  of  liquid  rocket  engine  heat  transfer  which  involve  the  injector  faceplate 
and  regeneratively  and  film  cooled  walls  are  being  investigated  by  computational  analysis.  A 
conjugate  heat  transfer  analysis  will  be  used  to  describe  localized  heating  phenomena 
associated  with  particular  injector  configurations  and  coolant  channels  and  film  coolant 
dumps.  These  components  are  being  analyzed,  and  the  analyses  verified  with  appropriate  test 
data.  Finally,  the  component  analyses  will  be  synthesized  into  an  overall  flowfield/heat 
transfer  model.  The  FDNS  code  is  being  used  to  make  the  component  analyses.  Particular 
attention  is  being  given  to  the  representation  of  the  thermodynamic  properties  of  the  fluid 
streams  and  to  the  method  of  combining  the  detailed  models  to  represent  overall  heating.  Unit 
flow  models  of  specific  coaxial  injector  elements  have  been  developed  and  will  be  described. 
Film  cooling  simulations  of  film  coolant  flows  typical  of  the  subscale  STME  being 
experimentally  studied  by  Pratt  and  Whitney  have  been  made,  and  these  results  will  be 
presented.  Other  film  coolant  experiments  have  also  been  simulated  to  verify  the  CFD  heat 
transfer  model  being  developed  by  SECA.  The  status  of  this  entire  study  will  be  presented, 
and  its  relevance  as  a new  design  tool  will  be  demonstrated. 


SECA,  Inc.,  3313  Bob  Wallace  Avenue,  Suite  202,  Huntsville,  AL 
Engineering  Sciences,  Inc.,  4920  Corporate  Drive,  Suite  K,  Huntsville,  AL 
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k-EQ.  CORRECTION: 
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1.44(1  +0.08  M°  25 ),  Cc2  = 1.92 
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• HOLDEN’S  TEST  CASE  #45 
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Mesh:  101  x 81 
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APPLICATION  OF  COMPUTATIONAL  FLUID  DYNAMICS 
TO  THE  DESIGN  OF  THE  FILM  COOLED  STME  SUBSCALE  NOZZLE 
FOR  THE  NATIONAL  LAUNCH  SYSTEM. 


By 

Joseph  L.  Garrett 

Government  Engine  and  Space  Propulsion 
Pratt  and  Whitney 
West  Palm  Beach,  Florida 


ABSTRACT 

A status  of  CFD  calculations  for  the  STME  film/dump  cooled  nozzle  design  will  be 
presented,  with  an  emphasis  on  the  timely  impact  of  CFD  in  the  design  of  the  sub-scale  nozzle 
coolant  system.  The  following  aspects  of  the  sub- scale  coolant  delivery  system  were  analyzed 
with  CFD: 

1.  Design  trade  study  of  a mechanical  flow  splitting  device  for  uniform  distribution  of  the 
subsonic  cavity  flow. 

2.  Design  trade  study  of  the  subsonic  cavity  lip  to  achieve  coolant  film  integrity. 

3.  Analysis  of  the  primary  flow  interaction  with  the  core/secondary  coolant  streams. 

All  design  calculations  were  performed  with  the  Generalized  Aerodynamic  Simulation  Program 
(GASP),  a 3-D,  multi-block,  generalized  Navier-Stokes  code  capable  of  solving  with  frozen, 
finite-rate  or  equilibrium  chemical  kinetics. 

The  initial  design  of  the  subsonic  cavity  flow  used  square  posts  to  distribute  the  sonic  orifice 
jets  into  a uniform  flow.  Calculations  for  this  design  indicated  that  an  unacceptable 
mal-distribution  of  film  occurred.  Design  modifications  involving  curved  and  slotted  posts  were 
computed  in  an  effort  to  uniformly  distribute  the  secondary  coolant  flow.  Analysis  of  these 
configurations  showed  that  although  the  flowfield  improved  in  uniformity,  it  was  still 
unacceptable,  especially  at  higher  feed  pressures.  Results  from  these  studies  were  then 
incorporated  into  a design  that  resulted  in  the  insertion  of  a porous  metal  ring  into  the  subsonic 
cavity.  Subsequent  water  flow  model  studies  showed  that  this  concept  was  successful  in 
uniformly  distributing  flow  exiting  the  cavity. 

In  addition  to  the  design  of  the  subsonic  cavity,  CFD  was  also  used  to  analyze  the 
secondary  coolant  lip  and  the  primary  flow  interaction  with  the  core/secondary  coolant  streams. 
A series  of  calculations  were  first  performed  to  modify  the  subsonic  cavity  lip  contour.  The  flow 
over  the  modified  lip  was  then  computed  simultaneously  with  the  primary  injectors  to  determine 
the  impact  of  the  subsonic  coolant  stream  on  the  primary  slot  jets. 

Pressure,  temperature,  velocity  and  coolant  mass  fraction  contours  will  be  presented  for 
these  configurations. 
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ANALYSIS  OF  SECONDARY  COOLANT  CAVITY 

Curve  Post  Design,  Low  P 
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ANALYSIS  OF  SECONDARY  COOLANT  CAVITY 

Slotted  Post  Design,  Low  P 
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ANALYSIS  OF  SECONDARY  COOLANT  CAVITY 

Comparison  of  Exit  Solutions 
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ANALYSIS  OF  SECONDARY  COOLANT  CAVITY 

Design  Solution:  Porous  Stainless  Steel  Filter 
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SUBSCALE  CORE/FILM  COOLANT  INTERACTION 

Previous  Assumptions 
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CORE/SECONDARY  COOLANT  FLOW  INTERACTION 

Velocity  Vectors  / H2  Mass  Fraction  Contours  Show  Vortex  Mixing 
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COMPARISON  OF  OLD  & NEW  SECONDARY  SLOT 

Pertinent  Geometry  And  Inflow  Conditions 
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COMPARISON  OF  OLD  & NEW  SUBSONIC  SLOT 

New  Conditions  Result  In  Substantial  Improvement  in  Film  Effectiveness 
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Note:  Increasing  Lip  Radius  and  Flow  Rate  Result  in  Diminishing  the  Impact  of  Vortex  Mixing 


COMPARISON  OF  OLD  & NEW  SUBSONIC  SLOT 

New  Conditions  Indicate  Dramatic  Improvement  In  Film  Cooling 
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Note:  Backside  Cooling  Over  The  Last  1/3  of  Injector  Ring  Would  Further  Reduce  Temperature  Load 


SUBSCALE  CORE/FILM  COOLANT  INTERACTION 

Imposed  Boundary  Values  For  High  m Case 
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MESH  SIZE:  75x61,  25x69,  59x81,  36x118,  36x118;  19,575  Grid  Points 


SUBSCALE  CORE/FILM  COOLANT  INTERACTION 

New  Conditions  Produce  Less  Mixing  Of  Secondary  Film 
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SUBSCALE  CORE/FILM  COOLANT  INTERACTION 

Secondary  Film  Layer  Provides  Adequate  Cooling  for  Injector  Ring 
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SUBSCALE  CORE/FILM  COOLANT  INTERACTION 

Integrity  of  Primary  Jet  Maintained 


Note:  Mixing  of  Jet  with  the  core  Indicates  Exhaust  Film  is  Less  than  50%  H2 


SUBSCALE  CORE/FILM  COOLANT  INTERACTION 

Nozzle  Geometry  & Film  Injection  Produce  Shock  Losses 
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Note:  Performance  Losses  to  be  Quantified 
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N92-3?255 


COMPUTATIONAL  FLUID  DYNAMICS  ANALYSIS  OF  SPACE 
SHUTTLE  MAIN  ENGINE  MULTIPLE  PLUME  FLOWS  AT 
HIGH-ALTITUDE  FLIGHT  CONDITIONS 

by 


N.  S.  Dougherty,  J.  B.  Holt,  B.  L.  Liu  and  S.  L.  Johnson 
Rockwell  International 
Space  System  Division 
Huntsville,  Alabama  35806 

ABSTRACT 


for  the  Consideration  of 

WORKSHOP  FOR  COMPUTATIONAL  FLUID  DYNAMIC 
APPLICATIONS  IN  ROCKET  PROPULSION 
APRIL  28-30, 1992 


Computational  fluid  dynamics  (CFD)  analysis  is  providing  verification  of 
Space  Shuttle  flight  performance  details  and  is  being  applied  to  Space 
Shuttle  main  engine  multiple  plume  interaction  flow  field  definition. 
Advancements  in  real-gas  CFD  methodology  described  herein  have  allowed 
definition  of  exhaust  plume  flow  details  at  Mach  3.5  and  107,000  ft.  The 
specific  objective  of  the  study  includes  the  estimate  of  flow  properties  at 
oblique  shocks  between  plumes  and  plume  recirculation  into  the  Orbiter 
base  so  that  base  heating  and  base  pressure  can  be  model  accurately.  The 
approach  utilizes  the  Rockwell  USA  Real  Gas  Three-Dimensional  Navier- 
Stokes  (USARG3D)  Code  for  the  analysis.  The  code  has  multi-zonal 
capability  to  detail  the  geometry  of  the  plumes  and  base  region  and  utilizes 
finite-rate  chemistry  to  compute  the  plume  expansion  angle  and  relevant 
flow  properties  at  altitude  correctly.  Through  an  improved  definition  of  the 
base  recirculation  flow  properties,  heating  and  aerodynamic  design 
environments  of  the  Space  Shuttle  Vehicle  can  be  further  updated. 

Results  of  IRAD  work  in  progress  indicate  that  at  this  altitude  the  plumes 
intersect  and  produce  oblique  shocks.  At  the  hottest  spot  of  the  oblique 
shock  between  Engines  2 and  3,  the  recovery  pressure  and  temperature 
were  found  to  be  216  psfa  and  5000  °F.  There  the  flow  resembles  a location  of 
a source  flow  with  a discriminating  streamline  driving  hot  gas  back  into 
the  base.  Considerable  exhaust  gas  is  forced  back  toward  the  thermal 
shield  of  the  Orbiter  base  between  Engines  2 and  3,  although  Engine  1 flow 
is  aspirated.  The  highest  base  temperature  at  the  lower  center  of  the  heat 
shield  reaches  2500  °F.  Future  work  is  planned  to  integrate  the  base  flow 
solution  to  the  integrated  vehicle  forebody  flow  for  a total  'nose-to-plume' 
solution.  With  the  vertical  tail  and  OMS  pods  effects  included  later,  it  is 
expected  that  there  may  be  recirculation  of  Engine  1 flow,  also. 
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Computational  Fluid  Dynamics  Analysis  of 
Space  Shuttle  Vehicle  and  Exhaust  Plume  Flows 
at  High  Altitude  Flight  Conditions 
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Rockwell  SPACE  SHUTTLE  SURFACE  GRID 

International 

Space  Systems  Division  Huntsville  Operations 
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N92-32256 


Direct  Numerical  Simulation  of  a Combusting  Droplet 

with  Convection 

P.  Y.  Liang 

CFD  Technology  Center, 

Rocketdyne  Division,  Rockwell  International 
Canoga  Park,  California 

The  evaporation  and  combustion  of  a single  droplet  under  forced  and  natural 
convection  has  been  studied  numerically  from  first  principles  using  a numerical 
scheme  that  solves  the  time-dependent  multi-phase  and  multi-species  Navier-Stokes 
equations  and  tracks  the  sharp  gas-liquid  interface  cutting  across  an  arbitrary  Eulerian 
grid.  The  flow  fields  both  inside  and  outside  of  the  droplet  are  resolved  in  a unified 
fashion.  Additional  governing  equations  model  the  inter-phase  mass,  energy  and 
momentum  exchange.  Test  cases  involving  iso-octane,  n-hexane  and  n-propanol 
droplets  show  reasonable  comparision  with  experimental  data  regarding  parameters 
such  as  breakup  mode,  evaporation  rate  and  flame  stand-off  distance.  The  partially 
validated  code  is  thus  readied  to  be  applied  to  more  demanding  droplet  combustion 
situations  where  substantial  drop  deformation  render  classical  models  inadequate. 
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DIRECT  NUMERICAL  SIMULATION  OF  A 


Control  Volume  Next  to  Droplet  Surface 


Rocketdyne  Division 
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Governing  Equations 
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ARICC-ST  MODELS 
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Fig.  1 Schematic  of  Partitioned  Computational  Volume  with 
Evaporating  Free  Surface  Segments 


Mode  Breakup  of  Simulated  (a  & c) 
Experimental  (b)  Droplet  Due  to 
k Wave  Passage,  (c)  Mas  done  with 
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Table  1.  Highlight  of  Problem  Parameters  for  Combusting  Droplet 
With  Convection  Simulations 
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t (sec) 

Fig.  3 Numerical  vs  Experimental  Drop  Diameter 
and  Drop  Center  Temperature  Histories 
in  Forced  Convection 
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Sequence  of  Simulated  Tempera lure  Contours  in  Grey  Scale 
Showing  Initial  Transient  of  n-Propanol  Droplet  Combusting 
with  Natural  Convection,  (white  Tmax=2500  K,  black  Tmin=30U  K) 
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Fig. 8 Close-up  View  of  Velocity  Vector  Plot  from  n-Propanol 
Drop  Combustion  Under  Natural  Convection  Simulation 
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Log^  Species  Concentration  Contours  from  n-Propanol  Drop  Combustion 
Under  Natural  Convection  Simulation 


CONCLUSIONS 
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LARGE  GRID  CONCENTRATION  AROUND  BOUNDARY  LAYER 
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A NUMERICAL  MODEL  FOR  ATOMIZATION-SPRAY 
COUPLING  IN  LIQUID  ROCKET  THRUST  CHAMBERS 

ty 

M.  G.  Giridharan,  A.  Krishnan,  J.  J.  Lee  and  A.  J.  Przekwas 
CFD  Research  Corporation,  Hunstville,  AL  35802 

and 

K.  Gross 

NASA  Marshall  Space  Flight  Center,  Hunstville,  AL  35802 


The  physical  process  of  atomization  is  an  important  consideration  in  the  stable 
operation  of  liquid  rocket  engines.  Many  spray  combustion  CFD  codes  do  not 
include  an  atomization  sub-model  but  assume  arbitrary  drop  size  distributions,  drop 
initial  locations  and  velocities.  It  has  been  shown  that  the  results  of  spray 
combustion  models  are  extremely  sensitive  to  the  assumed  droplet  initial 
conditions.  Furthurmore,  the  atomization  process  itself  is  a strong  function  of  the 
local  conditions  of  the  liquid  and  gas  flow.  Thus  it  is  important  to  account  for  the 
strong  mutual  coupling  between  the  liquid  phase,  the  spray  dynamics  and  the  gas 
flow.  A method  of  coupling  an  atomization  model  with  the  spray  model  in  a 
REFLEQS  CFD  code  will  be  presented.  This  method  is  based  on  a novel  Jet- 
Embedding  technique  in  which  the  equations  governing  the  liquid  jet  core  are 
solved  separately  using  the  surrounding  gas  phase  conditions.  The  droplet  initial 
conditions  are  calculated  using  a stability  analysis  appropriate  for  the  atomization 
regime  of  liquid  jet  break-up. 

This  novel  coupling  model  is  used  to  analyze  the  SSME  fuel  preburner  single 
injector  flow.  Results  of  the  diffusion  flame  characteristics  in  a single  injection 
element  will  be  presented.  The  effect  of  relative  velocity,  mixture  ratio  and  droplet 
initial  conditions  will  be  shown.  The  predictions  of  present  atomization  model  is 
compared  with  that  of  the  widely  used  CICM  correlation.  The  results  are  also 
compared  with  the  predictions  of  volume-of-fluid  method. 
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Spray  Model 

- Eulerian/Lagrangian  Particle  Tracking 
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Oxygen  Mass  Source  (gms/sec)  Distribution 


Ug  = 200.0  m/s,  Ul  = 25.0  m/s,  Tg  = 200  K 
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Oxygen  Mass  Source  ( gms/sec)  Distribution 
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Passive  Scalar  ( Hydrogen)  Concentration  Distribution 
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Numerical  Modeling  for  Dilute  and  Dense  Sprays 


C.P.  Chen.  Y.M.  Kim,  H.M.  Shang,  and  J.P.  Ziebarth 
University  of  Alabama  in  Huntsville 

T.S.  Wang 

NASA  Marshall  Space  Flight  Center 


Abstract 

Numerical  modelings  of  fuel-droplet  spray  combustion  finds  useful  applications 
for  the  assessment  of  the  engine  performance  stability  characteristics.  With  our  on- 
going studies  on  turbulent  reacting  flows,  we  have  successfully  implemented  a numer- 
ical model  for  spray-combustion  calculations.  In  this  model,  the  governing  gas-phase 
equations  in  Eulerian  coordinate  are  solved  by  a time-marching  multiple  pressure 
correction  procedure  based  on  the  operator-splitting  technique.  The  droplet-phase 
equations  in  Lagrangian  coordinate  are  solved  by  a stochastic  discrete  particle  tech- 
nique. In  order  to  simplify  the  calculation  procedure  for  the  circulating  droplets,  the 
effective  conductivity  model  is  utilized.  This  vaporization  model  includes  the  effects 
of  variable  thermophysical  properties,  non-unitary  Lewis  number  in  the  gas-film,  the 
Stefan  flow  effect,  and  the  effect  of  internal  circulation  and  transient  liquid  heating. 
The  k — e models  are  utilized  to  characterize  the  time  and  length  scales  of  the  gas- 
phase  in  conjunction  with  turbulent  modulation  by  droplets  and  droplet  dispersion 
by  turbulence.  This  method  entails  random  sampling  of  instantaneous  gas  flow  prop- 
erties and  the  stochastic  process  requires  a large  number  of  computational  parcels  to 
produce  the  satisfactory  dispersion  distributions  even  for  rather  dilute  sprays. 

The  present  study  has  made  two  major  improvements  in  spray  combustion  mod- 
elings. Firstly,  we  have  developed  a probability  density  function  approach  in  multi- 
dimensional space  to  represent  a specific  computational  particle.  Advantages  of  a 
parcel  PDF  tracking  method  is  to  reduce  the  number  of  computational  parcels  rep- 
resenting the  spray  dynamics  as  well  as  to  obtain  grid-  independent  solutions  for 
two-phase  flows.  Secondly,  we  incorporate  the  TAB  Taylor  Analogy  Breakup  model 
for  handling  the  dense  spray  effects.  These  breakup  models  is  based  on  the  reasonable 
assumption  that  atomization  and  drop  breakup  are  indistinguishable  processes  within 
a dense  spray  near  the  nozzle  exit.  Accordingly,  atomization  is  prescribed  by  injecting 
drops  which  have  a characteristic  size  equal  to  the  nozzle  exit  diameter. 

Example  problems  include  the  nearly  homogeneous  and  inhomogeneous  turbulent 
particle  dispersion,  and  the  non-evaporating,  evaporating,  and  burning  dense  sprays. 
Comparison  with  experimental  data  will  be  discussed  in  detail. 
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combustion. 
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- Multiple  pressure  correction 

- Non-iterative  for  transient  calculation 

- Applicable  to  all-speed  flows 
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GROUP/WIDTH  DISPERSION  MODEL 
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Time  step  to  interact  with  kth  eddy: 
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existing  variance  ot  the  particle  p 


Cumulative  pdf  distribution  at  any  point  in 
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Numerical  Implementations  of  Collision  Model 

Mean  collision  rate  between  each  pair  of  ’’parcels”  coexisting  i 
given  numerical  cell: 
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random  velocity  components  normal  to  the  relative  velocity 
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Figure  1.  Particle  dispersion  of  a nearly-honiogeneous  flow  for 
SSF  model(5000  particles)  and  DDWT  model. 
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Figure  4.22  Spray  parcel  distribution  in  a solid-cone  spray(t  = 3.0 ms) 
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Figure  4.23  Spray  tip  penetration  versus  time  in  a solid-cone  spray 
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SAUTER  MEAN  DIAMETER 
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Figure  4.25  Spray  parcel  distribution  and  velocity  vectors  in  a hollow-cone  spray 
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Figure  4.2G  Spray  tip  penetration  versus  time  in  a hollow-cone  spray 
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Figure  4.27  Spray  parcel  distribution  and  contours  of  fuel  mass  fraction  in  an 
evaporating  spray 
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Figure  4.28  Spray  tip  penetration  versus  time  in  an  evaporating  spray 
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Figure  4.29  Spray  parcel  distribution  and  contours  of  fuel  mass  fraction  in  a burn- 
ing spray 
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N92- 32259 

Modeling  of  SSME  Fuel  Preburner  ASI 

P.  Y.  Liang 

CFD  Technology  Center 

Rocketdyne  Division,  Rockwell  International  Corporation 
Canoga  Park,  California 

The  Augmented  Spark  Ignitor  (ASI)  is  a LOX/H2/electrical  spark  system  that  functions 
as  an  ignition  source  and  sustainer  for  stable  combustion.  It  is  used  in  the  SSME  pre- 
burner combustor,  the  SSME  main  combustion  chamber,  the  J-1  and  J-2  engines  as 
well  as  proposed  designs  of  the  Space  Transportation  Main  Engine  (STME)  main 
combustor  and  gas  generators.  In  the  SSME  it  is  a long  circular  cylindrical  chamber 
located  along  the  Main  Combustor  centerline  with  a truncated  conical  dome  at  the  top, 
which  contains  two  oblique  LOX  injection  ports  and  two  spark  plugs  offset  at 
90  degrees.  Hydrogen  injection  is  through  a number  of  nearly  tangential  slots  down- 
stream which  creates  a swirl  flow  intended  to  cool  the  ASI  chamber  walls.  Past  inci- 
dents of  erosion  of  the  ASI  spark  plugs  have  often  led  to  the  need  for  replacement  of 
these  very  costly  devices.  Thus  it  is  desirable  to  understand  the  complex  reactive  flow 
field  within  the  ASI  both  during  the  initial  ignition  transient  and  during  the  main  stage 
steady  state  combustion  (no  sparking). 

While  it  is  impossible  to  perform  direct  optical  diagnostics  to  measure  the  internal  flow 
field  of  the  ASI  under  hot-fire  conditions,  recent  advances  in  CFD-based  combustion 
modeling  have  made  it  feasible  to  characterize  the  flow  through  time-accurate  simula- 
tions. This  paper  documents  an  undertaking.to  characterize  the  flow  of  the  ASI.  The 
code  consists  of  a marriage  of  the  Implicit-Continuous-Eulerian/Arbitrary-Lagrangian- 
Eulerian  (ICE-ALE)  Navier-Stokes  solver  with  the  Volume-of-Fluid  (VOF)  methodology 
for  tracking  of  two  immiscible  fluids  with  sharp  discontinuities.  Spray  droplets  are 
represented  by  discrete  numerical  parcels  tracked  in  a Lagrangian  fashion. 
Numerous  physical  sub-models  are  also  incorporated  to  describe  the  processes  of 
atomization,  droplet  collision,  droplet  breakup,  evaporation,  and  droplet  and  gas 
phase  turbulence.  An  equilibrium  chemistry  model  accounting  for  8 active  gaseous 
species  is  also  used.  Taking  advantage  of  the  symmetry  plane,  half  of  the  actual  ASI 
is  modeled  with  a 3-dimensional  grid  that  geometrically  resolves  the  LOX  ports,  the 
spark  plug  locations,  and  the  hydrogen  injection  slots.  The  pertinent  features  and  for- 
mulations of  these  submodels  will  be  briefly  described  in  the  paper. 


1013 


MODELING  OF  SSME  FUEL  PREBURNER  ASI 


Rockwell  International 

Rocketdyne  Division 


SIGNIFICANCE 


UJ 

5 

UJ 

o 


DO 


UJ 

> UJ 

go 

o< 

u-  CO 

H 

< Q 
Ul  2 
DC< 

Q h- 
Z Z 

< UJ 

i 


m — 

<o 

DC  Z 

CO  DC 
UJ  3 
Q Q 


Q-  h- 
^ Z 

3* 


fiUJ 
O dc 
UJ  3 
O CO 

<< 

— J UJ 

o.  2 

Ml 

£ uj 


Rockwell  International 

Rocketdyne  Division 


(Dll 


1019 


Comparision  of  computed  vs  experimental 
drop  penetration  distances  for  different 
transverse  gas  velocities. 


orifice  diameter  d. 

Jo 


Fig.  2 Variation  of  penetration  with  orifice 
diameter;  Vg= 1 00  m/s,  V ^ = 1 5 m/s. 
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Fig.  1 Schematic  of  SSME  Preburner  ASI.  (LOX  ports  are 

actually  90  deg.  apart  from  spark  plugs  rather  than 
being  on  same  meridian  plane  as  shown) 


1023 


2 Close-up  Side  and  Top  Views  of  ASI  Dome  Region 


k = 1 — 


k = 1 1 — 
k = 1 5 — 


Fig.  3 Computational  Grid  (Partial)  Across  a J-plane  (shown 

reflected  across  central  axis  to  represent  full  chamber) 
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(a)  j =7  plane  (b)  k = 6 plane 

Fig.  4 Velocity  Profiles  In  Vertical  and  Horizontal  Cross 
Sectional  Planes  Before  Ignition  (t=.737  msec) 
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(b)  velocity  at  j=7 


(c)  temperature  at  k=6 

Fig.  5 Transient  Temperature  and  Velocity  Fields  During  and  Right 
After  Ignition  by  Spark  Plugs 
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(a)  velocity  (b)  temperature  (c)  log  0^  cone. 

Fig.  6 Velocity,  Temperature  and  Log  0^  Concentration  Profiles 
Across  j=7  Plane  at  Pseudo-steady  State  (t=4.10  msec) 


(a)  velocity 

Fig.  7 Velocity, 
Across  k= 


(b)  temperature  (c)  log  0,,  cone. 


Temperature  and  Log  02*  Concentration  Profiles 
6 Plane  at  Pseudo-steaay  State  (t=4.10  msec) 
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SUMMARY  OF  OBSERVATIONS  ON  ASI 
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CFD  MODELING  OF  TURBULENT  FLOWS  AROUND  THE  SSME  MAIN 
INJECTOR  ASSEMBLY  USING  POROSITY  FORMULATION 


Gary  C.  Cheng*,  Y.S.  Chen\  and  Joseph  H.  Rufi 


Abstract 


Hot  gas  turbulent  flow  distributions  around  the  main  injector  assembly  of  the  Space 
Shuttle  Main  Engine  (SSME)  and  LOX  flow  distributions  through  the  LOX  posts  have  great 
effect  on  the  combustion  phenomenon  inside  the  main  combustion  chamber.  An  advanced 
computational  fluid  dynamics  (CFD)  analysis  will  help  to  provide  more  accurate  and  efficient 
characterization  of  this  type  of  flow  field.  In  order  to  design  a CFD  model  to  be  an  effective 
engineering  analysis  tool  with  good  computational  turn-around  time  (especially  for  3-D  flow 
problems)  and  still  maintain  good  accuracy  in  describing  the  flow  features,  the  concept  of 
porosity  is  employed  to  describe  the  effects  of  blockage  and  drag  force  due  to  the  presence  of 
the  LOX  posts  in  the  turbulent  flow  field  around  the  main  injector  assembly  of  the  SSME.  A 
validated  non -isotropic  porosity  model  is  developed  and  incorporated  into  an  existing  Navier- 
Stokes  flow  solver  (FDNS).  Volume  and  surface  porosity  parameters,  which  are  based  on  the 
configurations  of  local  LOX  post  clustering,  will  be  introduced  in  to  the  governing  equations, 
which  can  be  written  as 


]_ 

J 


3pq 

9t 


3v.F. 

> l 

■sir 


+ V 


So  + R 

<,  q 


where  J is  the  Jacobian,  F;  is  the  sum  of  the  convective  flux  and  the  viscous  flux,  v;  is  the 
surface  porosity,  vv  is  the  volume  porosity,  Sq  and  Rq  are  the  source  and  residual  terms  of  the 
flow  variable  q,  respectively.  The  drag  force  and  the  heat  flux  source  due  to  the  presence  of 
LOX  posts  will  be  added  to  the  residual  term.  2-D  numerical  studies  have  been  conducted  to 
identify  the  drag  coefficients  of  the  flows  both  through  tube  banks  and  around  the  shielded  posts 
with  a wide  range  of  Reynolds  numbers.  A verified  model  of  the  drag  coefficients  is 
incorporated  into  the  FDNS  flow  solver.  A 2-D  flow  study  of  the  main  injector  assembly  is 
performed  to  verify  the  proposed  porosity  model.  A reasonable  O/F  ratio  distribution  was 
obtained,  therefore,  a 3-D  CFD  analysis  is  conducted  with  confidence.  The  3-D  CFD  analysis 
of  the  SSME  main  injector  assembly  is  divided  into  three  parts,  LOX  dome,  LOX  post  assembly 
torus,  and  hydrogen  cavity.  A62x91  x 16  mesh  system  is  constructed  for  the  LOX  dome, 
where  a 37  x 91  x 25  grid  system  is  employed  for  the  torus  region,  and  the  hydrogen  cavity  is 
discretized  into  a29x91  x 14  mesh  system.  The  numerical  study  of  the  turbulent  flow  in  the 
SSME  Phase  n+  power  head  is  analyzed  based  on  104%  power  balance  level,  and  the  result  will 
be  presented  in  the  coming  CFD  workshop  meeting. 


SECA,  Inc.,  3313  Bob  Wallace  Ave.,  Suite  202,  Huntsville,  AL 
Engineering  Sciences,  Inc.,  4920  Corporate  Dr.,  Suite  K,  Huntsville,  AL 
ED  32,  NASA/Marshall  Space  Flight  Center,  Huntsville,  AL 
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CFD  MODELING  OF  TURBULENT  FLOWS  WITHIN 
SSME  MAIN  INJECTOR  ASSEMBLY  USING  A 

POROSITY  MODEL 
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PRESENTATION  OVERVIEW 
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OBJECTIVE 
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3-D  POROSITY  MODEL/CFD  ANALYSIS  OF  THREE 
COMPONENTS  OF  THE  POWER  HEAD  (LOX  DOME,  LOX 
POST  ASSEMBLY,  HYDROGEN  CAVITY) 


PROPOSED  POROSITY  MODELS 
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SHIELDED  ELEMENTS 


a 

II 


o 


06 

ii 


O 


▲ ▲ 


o 


Q) 

I 

.O 


CL 


T 

CM 

O 

X 

• 

1 .785  x 1 04 

O 

X 

CD 

O 

00 

LO 

AP(psi) 

575.07 

625.49 

575.07 

rh  (Ib/sec) 

665.105 

105.65 

56.154 

Non-Baffle 

Elements 

Baffle  Elements 

First  Three  Rows 

1040 


LOX  POST  ASSEMBLY  POROSITY  MODEL  (104%  RPL) 
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3-D  POROSITY/CFD  ANALYSIS  OF  PHASE  11+  POWER  HEAD 


INCOMPRESSIBLE,  SINGLE  SPECIES 
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LOX  DOME , VELOCITY  VECTORS  AT  -30  DEG  (FUEL  SIDE)  PLANE 
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LOX  DOME , VELOCITY  VECTORS  AT  -2.  DEG 
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LOX  DOME , VELOCITY  VECTORS  AT  90  DEG  (OXIDIZER  SIDE)  PLANE 
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LOX  DOME,  VELOCITY  VECTORS  AT  178  DEG  PLANE 
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LOX  Done*  Exit  Pressure  Contours  (psf ) 


MAIN  INJECTOR  ASSEMBLY  WITH  TRANSFER  DUCTS 
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Geometry  of  Hot  Gas  Injector  Assembly 
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Fuel  Side 


Geometry  of  Hot  Gas  Injector  Assembly 
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ASSEMBLY  TORUS,  VELOCITY  VECTORS 
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ASSEMBLY  TORUS , VELOCITY  VECTORS  NEAR  BOTTOM  PLANE 


LOX  POST  ASSEMBLY  TORUS , VELOCITY  VECTORS  NEAR  MID  PLANE 
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LOX  Post  flssenbly.  Hot  Gas  Exit  Uelocity  Contours  (ft/sec) 
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LOX  Past  Rsseebly,  Hot  Gas  Exit  Pressure  Contours  (psf) 
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HYDROGEN  CAVITY 
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SYSTEM  OF  THE  HYDROGEN  CAVITY 
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Fuel  Side 


Hydrogen  Cai/ity,  Exit  Flou  Uelocity  Contours  (ft/sec) 
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0/F  Ratio  Contours  of  the  Main  Injector  Rssenbly  (1D4I  RPL ) 
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OF  THE  INJECTOR  FACE  (WITH  BLC  COOLANT  ADDED) 
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USING  PROPER  INLET  FLOW  PROFILES  TO  THE  LOX  DOME 


N92-322S1 

Computational  Fluid  Dynamics  Analysis 
of  SSME  Phase  II  and  Phase  11+  Preburner 
Injector  Element  Hydrogen  Flow  Paths 


Joseph  H.  Ruf 

Computational  Fluid  Dynamics  Branch 
Marshall  Space  Flight  Center,  NASA 


Phase  11+  Space  Shuttle  Main  Engine  powerheads  E0209  and 
E0215  degraded  their  Main  Combustion  Chamber  (MCC)  liners  at  a 
faster  rate  than  is  normal  for  phase  II  powerheads.  One  possible 
cause  of  the  accelerated  degradation  was  a reduction  of  coolant  flow 
through  the  MCC.  Hardware  changes  were  made  to  the  preburner 
fuel  leg  which  may  have  reduced  the  resistance  and,  therefore, 
pulled  some  of  the  hydrogen  from  the  MCC  coolant  leg. 

The  preburner  injector  element's  hydrogen  flow  path  changed 
significantly  from  the  phase  II  to  the  11+  design.  The  hydrogen  inlet 
area  was  reduced  by  42  percent,  the  annulus  length  was  shortened, 
and  the  geometry  of  the  annulus  convergence  section  changed.  With 
the  large  reduction  of  inlet  area,  an  increased  resistance  would 
normally  be  expected.  However,  a 10  percent  decrease  in  fuel  flow 
resistance  was  quoted  for  phase  11+  preburner  injector  elements. 

To  resolve  this  discrepancy,  a Computational  Fluid  Dynamics 
analysis  was  performed  to  determine  hydrogen  flow  path  resistances 
of  the  phase  11+  fuel  preburner  injector  elements  relative  to  the 
phase  II  element.  The  analysis  was  performed  for  104  percent  RPL 
conditions.  FDNS  was  implemented  on  axisymmetric  grids  with  the 
hydrogen  assumed  incompressible.  The  analysis  was  performed  in 
two  steps.  The  first  isolated  the  effect  of  the  different  inlet  areas 
and  the  second  modeled  the  entire  injector  element  hydrogen  flow 
path. 

The  isolated  effect  of  the  reduced  inlet  area  was  a 3.  percent 
increased  resistance  for  phase  11+  elements.  However,  the  entire 
flow  path  model  showed  no  difference  between  phase  II  and  11+ 
injector  element  resistances.  Phase  11+  annulus  geometry  changes 
compensated  for  the  reduced  inlet  area  such  that  there  was  no  net 
effect  on  hydrogen  flow  path  resistance. 
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CONCLUSION 


George  C.  Marshall  Space  Flight  Center 
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Computational  Fluid  Dynamics  Branch 


1074 


Analytically  determine  the  hydrogen  flow  path  resistance 
of  the  phase  11+  preburner  injector  elements. 


George  C.  Marshall  Space  Flight  Center 
Structures  and  Dynamics  Laboratory 
Computational  Fluid  Dynamics  Branch 
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has  a reduced  hydrogen  inlet  area  and  a shorter  hydrogen 
annulus. 

- Rocketdyne  quoted  a 1 0%  decrease  in  hydrogen  flow 
resistance  for  the  phase  11+  element  with  respect  to 
phase  II. 
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COOLANT  INLET 


PHASE  11+  INVESTIGATION 
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PREBURNER  INJECTOR  ELEMENTS 
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Assumptions: 
incompressible  flow 

rows  of  inlet  holes  modeled  as  equivalent  area 
circumferential  channels 
LOX  flow  not  modeled 
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Majority  of  the  flow  passes  through  the  holes  closest 
to  the  exit. 
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Row  Number 
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N92-32262 


STME  HYDROGEN  MIXER  STUDY 

Rob  Blumenthal 
Dongmoon  Kim 
George  Bache’ 


ABSTRACT 

The  hydrogen  mixer  for  the  STME  is  used  to  mix  cold  hydrogen  bypass  flow  with 
warm  hydrogen  coolant  chamber  gas,  which  is  then  fed  to  the  injectors.  It  is  very 
important  to  have  a uniform  fuel  temperature  at  the  injectors  in  order  to  minimize  mixture 
ratio  problems  due  to  the  fuel  density  variations.  In  addition,  the  fuel  at  the  injector  has 
certain  total  pressure  requirements.  In  order  to  achieve  these  objectives,  the  hydrogen 
mixer  must  provide  a thoroughly  mixed  fluid  with  a minimum  pressure  loss.  The 
AEROVISC  CFD  code  was  used  to  analyze  the  STME  hydrogen  mixer,  and  proved  to  be 
an  effective  tool  in  optimizing  the  mixer  design.  AEROVISC,  which  solves  die  Reynolds 
Stress-Averaged  Navier-Stokes  equations  in  primitive  variable  form,  was  used  to  assess 
the  effectiveness  of  different  mixer  designs.  Through  a parametric  study  of  mixer  design 
variables,  an  optimal  design  was  selected  which  minimized  mixed  rue!  temperature 
variation  and  fuel  mixer  pressure  toss.  The  use  of  CFD  in  the  design  process  of  the 
STME  hydrogen  mixer  was  effective  in  achieving  an  optimal  mixer  design  while  reducing 
the  amount  of  hardware  testing. 
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VIE  HYDROGEN  MIXER  PROVIDES  UNIFORM  TEMPERATURE 
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AN  EXPERIMENTAL  STUDY  OF  THE  FLUID  MECHANICS 
ASSOCIATED  WITH  POROUS  WALLS 
(AIAA  92-0769) 

N.  Ramachandran, 

Universities  Space  Research  Association, 

J.  Heaman, 

A.  Smith, 

NASA  Marshall  Space  Flight  Center, 

Huntsville,  A1  35812 

ABSTRACT 

The  fluid  mechanics  of  air  exiting  from  a porous  material  is  investigated. 
The  experiments  are  filter  rating  dependent,  as  porous  walls  with  filter 

ratings  differing  by  about  three  orders  of  magnitude  are  studied.  The  flow 
behavior  is  investigated  for  its  spatial  and  temporal  stability.  The  results  from 
the  investigation  are  related  to  jet  behavior  in  at  least  one  of  the  following 

categories:  (1)  Jet  coalescence  effects  with  increasing  flow  rate,  (2)  Jet  field 

decay  with  increasing  distance  from  the  porous  wall,  (3)  Jet  field  temporal 
turbulence  characteristics  and  (4)  Single  jet  turbulence  characteristics.  The 
measurements  show  that  coalescence  effects  cause  jet  development  and  this 
development  stage  can  be  traced  by  measuring  the  pseudoturbulence  (spatial 
velocity  variations)  at  any  flow  rate.  The  pseudoturbulence  variation  with 
increasing  mass  flow  reveals  an  initial  increasing  trend  followed  by  a 

leveling  trend,  both  of  which  are  directly  proportional  to  the  filter  rating.  A 
critical  velocity  begins  this  leveling  trend  and  represents  the  onset  of  fully 
developed  jetting  action  in  the  flow  field.  A correlation  is  developed  to  predict 
the  onset  of  fully  developed  jets  in  the  flow  emerging  from  a porous  wall.  The 
data  further  show,  that  the  fully  developed  jet  dimensions  are  independent  of 
the  filter  rating,  thus  providing  a length  scale  for  this  type  of  flow  field  (1 

mm).  Individual  jet  characteristics  provide  another  unifying  trend  with 
similar  velocity  decay  behavior  with  distance;  however,  the  respective 
turbulence  magnitudes  show  vast  differences  between  jets  from  the  same 
sample.  Measurement  of  the  flow  decay  with  distance  from  the  porous  wall 
show  that  the  higher  spatial  frequency  components  of  the  jet  field  dissipate 
faster  than  the  low  frequency  components.  Flow  turbulence  intensity 
measurements  show  an  out  of  phase  behavior  with  the  velocity  field  and  are 
generally  found  to  increase  as  the  distance  from  the  wall  is  increased. 
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"Experimental  Studies  of  Characteristic  Combustion-Driven  Flows  for  CFD  Validation" 

R.  J.  Santoro,  M.  Moser,  W.  Anderson,  S.  Pal,  H.  Ryan  and  C.  L.  Merkle 
Propulsion  Engineering  Research  Center 
The  Pennsylvania  State  University 
University  Park,  PA  16802 

A series  of  rocket-related  studies  intended  to  develop  a suitable  data  base  for  validation  of  CFD 
models  of  characteristic  combustion-driven  flows  has  been  undertaken  at  the  Propulsion 
Engineering  Research  Center  at  Penn  State.  Included  are  studies  of  coaxial  and  impinging  jet 
injectors  as  well  as  chamber  wall  heat  transfer  effects.  The  objective  of  these  studies  is  to 
provide  fundamental  understanding  and  benchmark  quality  data  for  phenomena  important  to 
rocket  combustion  under  well-characterized  conditions.  Diagnostic  techniques  utilized  in  these 
studies  emphasize  determinations  of  velocity,  temperature,  spray  and  droplet  characteristics  and 
combustion  zone  distribution.  Since  laser  diagnostic  approaches  are  favored,  the  development 
of  an  optically  accessible  rocket  chamber  has  been  a high  priority  in  the  initial  phase  of  the 
project.  During  the  design  phase  for  this  chamber,  the  advice  and  input  of  the  CFD  modeling 
community  were  actively  sought  through  presentations  and  written  surveys.  Based  on  this 
procedure,  a suitable  uni-element  rocket  chamber  has  been  fabricated  and  is  presently  under 
preliminary  testing.  Results  of  these  tests,  as  well  as  the  survey  findings  leading  to  the  chamber 
design,  was  presented. 

In  particular,  laser-induced  fluorescence  imaging  results  for  hydroxyl  radicals  have  been 
obtained.  These  experiments  were  conducted  using  gaseous  hydrogen/gaseous  oxygen 
propellants  in  the  optically  accessible  uni-element  rocket  chamber.  Heat  transfer  studies 
demonstrated  the  effectiveness  of  the  curtain  window  purge  needed  to  protect  the  optical 
surfaces.  These  results  also  demonstrated  the  capability  to  determine  wall  heat  transfer  rates  in 
the  present  rocket  test  chamber.  Measurements  indicated  that  heat  transfer  rates  were 
approximately  an  order  of  magnitude  smaller  than  observed  in  actual  rocket  chambers.  If  larger 
propellant  flow  rates  are  utilized,  this  difference  should  be  further  narrowed. 

Related  results  from  impinging  jet  studies  under  non-combusting  conditions  revealed  several 
trends  regarding  drop  size  and  liquid  sheet  breakup.  The  general  experimental  observation  that 
breakup  length  increases  with  increasing  jet  velocity  and  decreasing  impingement  angle  is  in 
agreement  with  previous  studies.  However,  this  trend  is  opposite  to  predictions  derived  from 
linear  stability-based  analysis  of  liquid  sheet  atomization.  However,  drop  size  predictions  for 
linear  stability -based  analysis  reproduced  the  observed  trend  of  decreasing  drop  size  with 
increasing  jet  velocity  and  increasing  impingement  angle. 
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OUTLINE 


• CFD  Validation  - Optically  Accessible  Rocket 

• Heat  Transfer  Measurements 


• Impinging  Jet  Studies 

• Summary 

• Future  Work 
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CFD  VALIDATION  - OPTICALLY 
ACCESSIBLE  ROCKET 

M.  Moser,  S.  Pal,  R.  Santoro,  C.  Merkle 


OBJECTIVE 


To  Provide  Benchmark-Quality  Data  for  CFD  Code 
Validation  for  Combustion-Driven  Flows 


Obtain  Fundamental  Data  Under  Realistic  and  Well 
Characterized  Conditions 
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APPROACH 


Emphasize  Uni-element  Coaxial  Injectors 


Obtain  Fundamental  Data  Under  Well  Characterized  and  Realistic 
Conditions 

Pressure 

Combustion  zone 
Mean  Velocity 
Turbulence  Intensity 
Species 

Droplet  Size  and  Velocity 


Employ  Non-Intrusive  Advanced  Diagnostics 

Laser  Induced  Fluorescence 

Raman  Spectroscopy 

Laser  Velocimetry 

Flow  Visualization 

Phase  Doppler  Particle  Analyzer 

Polarization  Ratio 

Raman  Spectroscopy 


Closely  Integrate  CFD  Objectives  Into  Experimental  Program 

Experiment  Design 

Measured  Quantities 

Boundary  Condition  Specification 
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SUMMARY  OF  SURVEY  RESULTS 


• Square  Geometry  Is  Acceptable 

Initially  Axisymmetric  Approximation  Suitable  Eventually 
3-d  Verification  Required 


• Chamber  Mach  Number  Need  Not  Be  Matched  to 
Actual  Rocket  Conditions 


• Recirculation  Effect  In  Uni-Element  Chamber  Can 
Be  Accommodated  If  Suitable  Measurements  Are 
Available 


• Measurements  And  Boundary  Conditions  - Nearly 
Everything  Is  Important 

Typically  All  Parameters  Rate  1 Or  2 On  Scale  Of  5 


Responses  From: 


MSFC 
Aeroj  et 
Rocketdyne 
CFDRC 
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OPTICALLY  ACCESSIBLE  ROCKET 


• Two  Inch  Square  Cross-section 

• Variable  Length  6-12  Inches 

• Two  Inch  Diameter  Quartz  Windows  For  Viewing 

• Slot  Windows  For  Laser  Access 

• Injector  Design  Is  Modular 
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RESULTS 

• Rocket  Chamber  Has  Been  Checked  Out  To  400 
psi  Chamber  Pressure 

• Window  Purges  Effectiveness  Tested  With  Heat 
Flux  Gages 


• Quartz  Windows  Tested  To  200  psi  For  Four 
Second  Firings 


• Shadowgraph  Photograph  Obtained 


• 2-D  Laser  Induced  Fluorescence  Image  Of  OH 

Obtained 
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as  sensors 
id  additional 
ita  systems 

Sequencer  and 
abort  monitor 


Control  panal 


Test  cell 


fiowcontrol 

panel 


Test  stand 
and  rocket 
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HEAT  TRANSFER  MEASUREMENTS 

S.  Pal,  C.  Merkle 

Objective: 

• Demonstrate  Capability 


• Test  Effectiveness  Of  Window  Purge 
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WINDOW  PURGE  COMPARISON 


,b 
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TIME  (s) 


OBJECTIVES 


O £ 
CO  CD 

CO  CD 

sz  c 
O § 

S''5 

2cl 

Q 

CO  « 

‘S.o 

>'’co 
• « 

> CD 
■•£=  Q. 

wO 

§T3 

W § 

3 O 
CO  'C= 

.0)q) 

w E 

CD  O 

> CD 

■^O 


*-  c c 

CD  <<  <D 

aC  E 
E c ® 

= CD  rr> 


1152 


• Compare  Experimental  Results  with  Theoretical 
Models  of  Sheet  Breakup  and  Drop  Formation 

- Linear  Growth  of  Surface  Waves 
Due  to  Aerodynamic  Forces 


Impinging  Jet  Injector  System 
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WATER  FLOW  ^ SPRAY  FAN 


EIGHT  ADJUSTMENT 


Uj  (m/s) 
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SUMMARY 
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• Pre-lmpingemenl  Length  Has  Measurable  Cliect 
on  Breakup  Length  and  Drop  Size 


SUMMARY 


• Successful  Firings  of  Optically  Accessible  Rocket 
Test  Chamber  Achieved 


• Shadowgraph  of  Combustion  Zone  Obtained 

• Preliminary  Two-Dimensional  OH  Imaging  Results 
Obtained 


• Wall  Heat  Transfer  Rate  Measurement  Capability 
Demonstrated 


• Fundamental  Studies  of  Impinging  Sprays  Are 
Elucidating  Basic  Atomization  Mechanisms 


1159 


FUTURE  WORK 


Semi-Quantitative  Measurement  of  Relative  OH 
Concentration 


Velocimetry  (initially  single  point  measurements, 
eventually  2-D  velocimetry  measurements) 


Raman  Spectroscopy 
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N92-32265 

TURBINE  DISK  CAVITY  AERODYNAMICS  AND  HEAT  TRANSFER 

B.  V.  Johnson 
W A.  Daniels 

United  Technologies  Research  Center 
Hast  Hartford.  (T  06108 

Experiments  were  conducted  to  define  the  nature  of  the  aerodynamics  and  heat  transfer 
for  the  flow  within  the  disk  cavities  and  blade  attachments  of  a large-scale  model,  simulating 
the  SSME  turbopump  drive  turbines.  These  experiments  of  the  aerodynamic  driving 
mechanisms  explored  the  following:  ( 1 ) flow  between  the  main  gas  path  and  the  disk  cavities, 
(2)  coolant  flow  injected  into  the  disk  cavities,  (3)  coolant  density,  (4)  leakage  flows  through 
the  seal  between  blades,  and  (5)  the  role  that  each  of  these  various  flows  has  in  determining 
the  adiabatic  recovery  temperature  at  all  of  the  critical  locations  within  the  cavities.  The 
model  and  the  test  apparatus  provide  close  geometrical  and  aerodynamic  simulation  of  all 
the  two-stage  cavity  flow  regions  for  the  SSME  High  Pressure  Euel  'Turbopump  and  the 
ability  to  simulate  the  sources  and  sinks  for  each  cavity  flow. 

Carbon  dioxide  was  used  as  a trace  gas  for  constant  density  experiments  or  as  the 
simulated  “heavy  gas”  coolant.  Gas  samples  were  withdrawn  at  selected  locations  on  the 
rotating  and  stationary  surfaces  in  the  fore  and  aft  cavity  and  the  interstage  seal  regions  of 
the  two  stage  system.  The  gas  samples  were  used  to  determine  the  fraction  of  gas  at  a location 
which  originates  from  each  of  three  coolant  injection  locations  or  four  gas  path  locations. 
Samples  were  also  withdrawn  at  selected  locations  in  the  blade  shank  regions. 

A parametric  series  of  experiments  was  conducted  with  constant  density  fluids  and  an 
exploratory  series  of  experiments  was  conducted  with  CO2  as  the  simulated  coolant. 
Experimental  results  showed  (1)  the  variation  of  coolant  distribution  on  the  cavity  and  disk 
surfaces  as  a function  of  coolant  flow  ratio,  (2)  the  effects  on  the  coolant  distribution  for 
changes  in  the  coolant  inlet  distributions,  and  (3)  increased  mixing  of  coolant  with  the 
ingested  gas  when  a heavy  gas  (density  ratio  equal  1.5)  was  used  as  the  coolant. 
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TURBINE  DISK  CAVITY  AERODYNAMICS 

AND  HEAT  TRANSFER 
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ACTUAL  AND  MODEL  DISK/CAVITY  SYSTEMS 
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MODEL  INSTRUMENTATION 
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□ Thermocouples  • Pressure/C02  taps  on  rotating  components 

O Pressure/C02  taps  in  passages  ▲ Pressure/CC>2  taps  on  stationary  components 


SSME  TURBINE  DISK  CAVITY  MODEL 


I 
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ORIGINAL  PAGE 

BLACK  AND  WHITE  PHOTOGRAPH 


C89-414-B-2 


MODEL  IN  INTERNAL  FLOW  FACILITY 


ORIGINAL 
BLACK  AND  Wh'if 


jtograph 


MODEL  SEAL  REGION  AND  GAS  SOURCE/EXIT  LOCATIONS 


COOLANT  CONCENTRATION  ON  ROTOR  AND  STATIONARY  WALLS 
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Dimensionless  coolant  flow  rate,  (mc/27T/iaR0)/(  Pa  OR02/  A#a)0-8 
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0.01  0.02  0 0.01  0.02 

Dimensionless  coolant  flow  rate,  (mc/2777/aR0)/(  pa  QR02/  /ia)0.8 


COOLANT  DISTRIBUTION  ON  ROTOR 
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COOLANT  DISTRIBUTION  ON  STATIONARY  WALL 
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Dimensionless  coolant  flow  rate,  ( mc/2777/aR0)/{  Pa  flR02/pa)0,8 


COOLANT  DISTRIBUTION  ON  ROTOR 


Dimensionless  coolant  flow  rate.  (mc/27T//aR0)/(  Pa  flR02/ /Ja)0*8 


COOLANT  DISTRIBUTION  ON  STATIONARY  WALL 
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RESULTS/CONCLUSIONS 
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RESULTS/SPECULATION 
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A Numerical  Study  of  Two-Dimensional  Vortex  Shedding 
From  Rectangular  Cylinders 

A.  H.  Hadid,  M.  M.  Sindir 
CFD  Technology  Center 
Rockwell  International/Rocketdyne  Division 
Canoga  Park.,  California 
R.  I.  Issa 

Department  of  Mineral  Resources  Engineering 
Imperial  College  of  Science 
Technology  and  Medicine, 

London,  SW7,  2BP,  England 

An  efficient  time-marching,  non  iterative  calculation  method  is  used  to  analyze  time- 
dependent  flows  around  rectangular  cylinders.  The  turbulent  flow  in  the  wake  region  of 
a square  section  cylinder  is  analyzed  using  an  anisotropic  k-e  model.  Initiation  and 
subsequent  development  of  the  vortex  shedding  phenomenon  is  naturally  captured 
once  a perturbation  is  introduced  in  the  flow.  Transient  calculations  using  standard 
eddy-viscosity  and  an  anisotropic  k-e  models  averaged  over  an  integral  number  of 
cycles  to  get  the  fluctuating  energy  (organized  and  turbulent)  are  compared  with 
experimental  data.  It  is  shown  that  the  anisotropic  k-e  model  resolves  the  anisotropy  of 
the  Reynolds  stresses  and  give  mean  energy  distribution  closer  to  the  experiment  than 
the  standard  k-e  model. 
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A NUMERICAL  STUDY  OF  TWO-DIMENSIONAL  VORTEX  SHEDDING  FROM 

RECTANGULAR  CYLINDERS 
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Rockwell  International 
Rocketdyne  Division 
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Rockwell  International 
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Rockwell  International 
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PISO"  NUMERICAL  PROCEDURE  (CONTD.) 
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MODEL  EQUATIONS 
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MODEL  EQUATIONS  (CONTD.) 


K 

* 2 
S § 

w o 
2 P 
w < 

CJ  ft- 

X M 

p 00 

E £2 

2 Q 

s >« 

z ^ 

a « 

i-J  2 

p W 

ffl  Q 
(X  w 
p o 

E-l  ■< 

Q & 

m W 

O > 
< < 
w 

W CO 

> < 

< 5 

U 

CO  W 

< P 

S H 

ft-  C/2 

w T. 


S2  =l 
£ u 
rl  11 

ftS  > 

S Q 
P 2 
£ < 


> l£ 

+ 

> 


I'D  * 

ro 


I'D  ns 


*-*  CO 
> & 


ns  * 
ns 


wi**- 

coins 


D Vf 
ns  ns 


1189 
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Cxm  (m=  1,2,3)  ARE  MODEL  CONSTANTS 


MODEL  EQUATIONS  (CONTD.) 
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Rockwell  International 
Rocketdyne  Division 


RESULTS  AND  DISCUSSIONS 
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RESULTS  AND  DISCUSSIONS  (CONTD.) 


1193 


Rockwell  International 
Rocketdyne  Division 


NORMAL  VELOCITY  HISTORY  AT  CENTERLINE  5 STEP  HEIGHTS 

DOWNSTREAM  AT  RE=14000 
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VELOCITY  ALONG  THE  CENTERLINE 


33333353*33^33 

on/n 


° -2 

§ -8  1 

’J  -a  i 

3 5 w 

•T'  **•  T 


5 3-1 
« " 1 5 
•JS  ,2  i s 

! 1 1 1 

-all 


1198 


Rockwell  International 

Rocketdyne  Division 


MEAN  KINETIC  ENERGY  ALONG  THE  CENTERLINE 


1199 


Rockwell  International 
Rocketdyne  Division 


MEAN  AXIAL  VELOCITY  ALONG  THE  CENTERLINE 


MEAN  KINETIC  ENERGY  ALONG  THE  CENTERLINE 
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NORMAL  STRESS  CONTOURS 


Rockwell  International 
Rocketdyne  Division 


RESULTS  OF  MEAN  KINETIC  ENERGY  SHOW  THAT  THE  ANISOTROPIC 
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for  the  CFD  Workshop  — 1992 
A Status  of  the  Turbine  Technology  Team  Activities 

Lisa  W.  Griffin 


A Status  of  the  activities  of  the  Turbine  Technology  Team  of  the 
Consortium  for  Computational  Fluid  Dynamics  (CFD)  Application  in 
Propulsion  Technology  is  presented.  The  team  consists  of 
members  from  the  government,  industry,  and  universities.  The 
goal  of  this  team  is  to  demonstrate  the  benefits  to  the  turbine 
design  process  attainable  through  the  application  of  CFD.  This 
goal  is  to  be  achieved  by  enhancing  and  validating  turbine 
design  tools  for  improved  loading  and  flowfield  definition  and 
loss  prediction,  and  transferring  the  advanced  technology  to  the 
turbine  design  process. 

In  order  to  demonstrate  the  advantages  of  using  CFD  early  in  the 
design  phase,  the  Space  Transportation  Main  Engine  (STME) 
turbines  for  the  National  Launch  System  (NLS)  were  chosen  on 
which  to  focus  the  team's  efforts.  The  Turbine  Team  activities 
run  parallel  to  the  STME  design  work. 

Work  during  the  past  year  has  centered  on  transferring 
technology  obtained  through  the  team ' s Generic  Gas  Generator 
Turbine  Program  (reported  on  in  the  1990  workshop)  to  the  STME 
LOX  turbine  design  point.  A preliminary  baseline  design  was 
analyzed  through  CFD,  and  areas  requiring  refinement  to 
eliminate  local  overspeeds  and  separations  were  found.  An 
improved  baseline  design  has  been  finalized.  The  team  is 
currently  comparing  results  from  five  team  members  (Aerojet, 
NASA/Ames  Research  Center,  NASA/ Lewis  Research  Center,  Pratt  & 
Whitney,  and  Scientific  Research  Associates).  These  solutions 
will  be  compared  to  data  to  be  obtained  in  the  Marshall  Space 
Flight  Center  Turbine  Airflow  Facility.  Interrogation  of  these 
solutions  are  also  in  progress  to  determine  high  loss  locations 
and  to  provide  guidance  for  developing  and/or  implementing 
concepts  to  control  these  losses. 
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New  Programs 

- Volute  Development 

- Supersonic  Turbine  Development 


A Summary  of  the  Activities  of  the  NASA/MSFC 
Turbine  Technology  Team 
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The  University  of  Alabama  (UA) 

The  University  of  Alabama  in  Huntsville  (UAH) 
Virginia  Polytechnic  institute  (VPI) 
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- Baseline  Unsteady  Aerodynamics/Heat  Transfer  (CUBRC) 

- Smooth  Blades  (MSFC) 

- Circumferential  Exit  AP  (MSFC) 

- ATD  (MSFC) 
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Focused  flow  analysis  tools  on  the  STME  LOX  turbine  design  point  resulting  in  Gas  Generator 
Oxidizer  T urbine  (GGOT) 
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Airfoil  Count  Base  -55  percent 


BASELINE  GGGT  AERODYNAMIC  DbSIGN 
Calculated  Streamlines  - Midspan 


CONSORTIUM  FOR  CFD  APPLICATION  IN  PROPULSION  TECHNOLOGY 

Turbine  Stage  Technology  Team-Baseline  GGGT  aerodynamic  analyses 

LeRC 


ORIGINAL  PAGE  IS 
OF  POOR  QUALITY 
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ANALYSIS  OF  A NEW  TURBINE  DESIGN 
Instantaneous  Temperature  Contours 
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- Development  of  concepts  to  control  losses 
- Parametric  studies  with  CFD 

- Development  of  advanced  design 

- Detailed  analyses  of  advanced  design 

- Experimental  evaluation  of  advanced  design 


CONSORTIUM  FOR  CFD  APPLICATION  IN  PROPULSION  TECHNOLOGY 
Turbine  Technology  Team  - Baseline  GGOT  Aerodynamic  Analyses 
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• SRA  - Steady  3D  NIS,Tip  Lkg 
•Aerojet  - Steady  3D  N/S,  Tip  Lkg 


GAS  GENERATOR  OXIDIZER  TURBINE  (GGOT)  DESIGN 

Streaklines  Downstream  of  Rotor 
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3D  Navier-Stokes 
Analyses, 

Steady  & Unsteady 


GGOT  ROTOR  - TIP  LEAKAGE  FLOW 
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GGOT  ROTOR  - PRESSURE  LOSS  COEFFICIENT  CONTOURS 
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(max  inlet)  - 


GGOT  ROTOR  - PRESSURE  LOSS  COEFFICIENT  CONTOURS 
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Experimental  evaluation  of  baseline  volutes  on  GGOT  tor  coae  vanaauon 
Development  of  loss  control  and  side  load  reduction  concepts 

— Parametric  studies  with  CFD 

~ Detailed  analysis  of  improved  designs 

- Experimental  evaluation  of  improved  volutes  on  advanced  concepts  GGOT 
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- Detailed  analyses  of  advanced  design 

- Experimental  evaluation  of  advanced  design 
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Team  activities  have  been  coordinated  and  focused  on  demonstration  of  CFD  tools  for  turbine 
design 


N92-32268 


A CRITICAL  EVALUATION  OF  A THREE-DIMENSIONAL 
NAVIER-STOKES  CFD  AS  A TOOL  TO  DESIGN 
SUPERSONIC  TURBINE  STAGES 


C.  HAH,  O.  KWON,  AND  M.  SHOEMAKER 
NASA  LEWIS  RESEARCH  CENTER 
21000  BROOKPARK  ROAD,  CLEVELAND,  OHIO 


Three-dimensional  flow  phenomena  in  a supersonic  turbine  blade  row  have  been 
studied  numerically  to  evaluate  CFD  as  a tool  to  design  supersonic  turbine  stages. 

The  details  of  the  three-dimensional  flow  structure  inside  the  supersonic  turbine 
blade  row  and  the  overall  aerodynamic  performance  at  design  and  off-design  condi- 
tions are  analyzed  and  the  results  are  compared  between  the  experimental  data  and 
the  numerical  results. 
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SUPERSONIC  TURBINE  STAGE 
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FLOWS  INSIDE  SUPERSONIC  TURBINE 


OBJECTIVES 
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FLOWS  INSIDE  SUPERSONIC  TURBINE 
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UNSTARTED  FLOW 


OVERALL  APPROACHES 
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COMPUTATIONAL  GRID 
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VELOCITY  VECTORS  FOR  UNSTARTED  FLOW 
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COMPARISON  FOR  STARTED  FLOW 
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MERIDIONAL  VELOCITY  VECTORS 


SHOCK-INDUCED  VORTEX  GENERATION 


CONCLUSIONS 
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N92-32269 

NAVIER-STOKES  ANALYSIS  OF  AN  OXIDIZER  TURBINE  BLADE  WITH  TIP  CLEARANCE* 


Howard  J.  Gibeling 
Jayant  S.  Sabnis* 
Scientific  Research  Associates,  Inc. 
Glastonbury,  CT 


ABSTRACT 

The  Gas  Generator  Oxidizer  Turbine  (GGOT)  Blade  is  being  analyzed  by  various 
investigators  under  the  NASA  MSFC  sponsored  Turbine  Stage  Technology  Team  design  effort. 
The  present  work  concentrates  on  the  tip  clearance  region  flow  and  associated  losses;  however, 
flow  details  for  the  passage  region  are  also  obtained  in  the  simulations.  The  present  calculations 
simulate  the  rotor  blade  row  in  a rotating  reference  frame  with  the  appropriate  coriolis  and 
centrifugal  acceleration  terms  included  in  the  momentum  equations.  The  upstream 
computational  boundary  is  located  about  one  axial  chord  from  the  blade  leading  edge.  The 
boundary  conditions  at  this  location  have  been  determined  by  Pratt  & Whitney  using  an  Euler 
analysis  without  the  vanes  to  obtain  approximately  the  same  flow  profiles  at  the  rotor  as  were 
obtained  with  the  Euler  stage  analysis  including  the  vanes.  Inflow  boundary  layer  profiles  are 
then  constructed  assuming  the  skin  friction  coefficient  at  both  the  hub  and  the  casing.  The 
downstream  computational  boundary  is  located  about  one  axial  chord  from  the  blade  trailing 
edge,  and  the  circumferentially  averaged  static  pressure  at  this  location  was  also  obtained  from 
the  P&W  Euler  analysis. 

Results  have  been  obtained  for  the  3-D  baseline  GGOT  geometry  at  the  full  scale  design 
Reynolds  number.  Details  of  the  clearance  region  flow  behavior  and  blade  pressure  distributions 
have  been  computed.  The  spanwise  variation  in  blade  loading  distributions  are  shown,  and 
circumferentially  averaged  spanwise  distributions  of  total  pressure,  total  temperature,  Mach 
number  and  flow  angle  are  shown  at  several  axial  stations.  The  spanwise  variation  of  relative 
total  pressure  loss  shows  a region  of  high  loss  in  the  region  near  the  casing.  Particle  traces  in 
the  near  tip  region  show  vortical  behavior  of  the  fluid  which  passes  through  the  clearance  region 
and  exits  at  the  downstream  edge  of  the  gap.  Future  work  will  include  collaboration  with  the 
P&W  design  team  to  identify  design  changes  which  may  reduce  clearance  flow  losses. 

* Work  supported  under  NASA  MSFC  Contract  NAS8-38865. 

* Currently  at  United  Technologies  Research  Center,  East  Hartford,  CT. 
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GGOT  SIMULATION  SUMMARY 
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25  POINTS  IN  CLEARANCE  REGION 
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OXIDIZER  TURBINE  BASELINE  DESIGN 

FULL  SCALE  TURBINE  FLOWPATH  CLOSE-UP 
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BASELINE  GGOT  HUB  SECTION  GRID 

(100x60) 
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BASELINE  GGOT  MID-SPAN  SECTION  GRID 

(100  x 60) 
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BASELINE  GGOT  TIP  SECTION  GRID 

(100  x 60) 
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BASELINE  GGOT  3-D  GRID  - HUB  AND  CASING  SECTIONS 
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BASELINE  GGOT  WITH  TIP  CLEARANCE 
MACH  NUMBER  CONTOURS  AT  9.2%  SPAN 
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BASELINE  GGOT  WITH  TIP  CLEARANCE 
MACH  NUMBER  CONTOURS  AT  90%  SPAN 
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BASELINE  GGOT  WITH  TIP  CLEARANCE 
PRESSURE  CONTOURS  AT  9.2%  SPAN 
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BASELINE  GGOT  WITH  TIP  CLEARANCE 
PRESSURE  CONTOURS  AT  53.5%  SPAN 
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Baseline  GGOT  with  Clearance 
Blade  Surface  Pressure  (9.2%  span) 
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Circumferentially  Averaged  Variables 
Total  Pressure  Loss  Across  Blade 
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PARTICLES  TRACES 
BASELINE  GGOT  WITH  TIP  CLEARANCE 
BLADE  TIP  STREAK  LINES 
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N92-32270 

NUMERICAL  SIMULATION  OF 
TURBOMACHINERY  FLOWS  WITH 
ADVANCED  TURBULENCE  MODELS 


B.  Lakshminarayana,  R.  Kunz,  J.  Luo,  S.  Fan 
The  Pennsylvania  State  University 


A three  dimensional  full  Navier-Stokes  (FNS)  code  is  used  to  simulate  complex 
turbomachinery  flows.  The  code  incorporates  an  explicit  multistep  scheme  and  solves  a 
conservative  form  of  the  density  averaged  continuity,  momentum  and  energy  equations  in 
body  fitted  coordinates.  Rotation  terms  are  included  for  the  computation  of  rotor  flows. 
A compressible  low-Reynolds  number  form  of  the  k-  e turbulence  model,  and  aq-u  model 
and  an  algebraic  Reynolds  stress  model  have  been  incorporated  in  a fully  coupled  manner 
to  approximate  Reynolds  stresses. 

The  code  is  used  to  predict  viscous  flow  field  in  a backswept  transonic  centrifugal 
compressor  for  which  laser  two  focus  data  is  available.  The  tip  clearance  flow,  and 
curvature  and  rotation  induced  secondary  flows  are  captured  to  good  accuracy.  Solutions 
which  incorporate  Reynolds  stress  model  show  significant,  though  not  dramatic,  differences 
in  predicted  secondary  flows,  wall  shear  stress  and  performance  parameters,  when  compared 
to  the  k-  e solution. 

The  code  is  also  utilized  to  simulate  the  tip  clearance  flow  in  a cascade.  An 
embedded  H-grid  topology  was  utilized  to  resolve  the  flow  physics  in  the  gap.  The  data  and 
predictions  were  performed  with  and  without  tip  clearance,  endwall,  wall  motion. 
Additionally,  both  Euler  and  Navier-Stokes  computations  were  performed.  The  results 
indicate  that  the  Navier-Stokes  code  captures  the  flow  physics  accurately,  including  tip 
vortex  strength,  trajectory,  loading  and  interaction  of  tip  clearance  flow  with  the  secondary 
flow. 


The  code  has  also  been  used  to  predict  the  two  dimensional  viscous  and  thermal  flow 
field  in  a transonic  turbine  nozzle  with  75°  turning,  = 015,  M2  = 0.7  to  1.11,  Re  = 0.5  - 
2 x 106,  and  T0  = 415°K.  Good  agreement  is  obtained  for  pressure  distribution,  wake  and 
surface  heat  transfer. 

The  code  has  been  extended  to  include  unsteady  Euler  solution  for  predicting  the 
unsteady  flow  through  a cascade  due  to  incoming  wakes,  simulating  rotor-stator  interaction. 
Unsteady  characteristic  boundary  conditions  are  specified  at  inlet  and  exit.  The  predicted 
unsteady  surface  pressure  distribution,  unsteady  lift  and  moment,  pressure  wave  propagation 
in  a flat  plate  due  to  an  incoming  gust  compares  well  with  the  analytical  theories  and  earlier 
computations.  This  code  will  be  integrated  with  a boundary  layer  code  to  capture  the 
unsteady  viscous  flow  and  heat  transfer. 


PRECEDING  PAGE  BLANK  NOT  FILMED 
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NUMERICAL  SIMULATION  OF 
TURBOMACHINERY  FLOWS  WITH 
ADVANCED  TURBULENCE  MODELS 


B.  Lakshminarayana,  R.  Kunz,  J.  Luo,  S.  Fan 
The  Pennsylvania  State  University 


1.  Introduction 

2.  Technique  and  Turbulence  Models 

3.  Computation  of  Transonic  Centrifugal  Compressor  Flow 
Field 

4.  Computation  of  Tip  Clearance  Flows  in  Cascades 

5.  Computation  of  Aerothermal  Field  in  a Transonic  Nozzle 

6.  Simulation  of  Rotor/Stator  Interaction 

7.  Conclusions 
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GOVERNING  EQUATIONS 
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Coakley's  q-co  model 


• q=Vk,  w=e/k 

* ^t=cnVpq2/“ 

cy=0.09 

f^=l-exp(-aR) 

a=0.02 


R= 


q-equation 

gjpq)  a(pqui)  a , , 3q 

at  + 9xj  axj[(,1+pt/Prq)ax  I 


+2(Vn  0)2  ' 3 CO 


‘J 

S 


2 V 


-1)pcoq 


3U:  3U:  dU: 

S= 

3Xj  3xj  5xj 

Prq=2.° 


p 3uj, 

3t> 


2 
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• co-equation 


3(pco)  ^(P®uj)  d dco 

at  + axj  =axjt(^+tJ-t/Pro))a7.] 

+[C1  (CH  ^ ' 3 P -c2lP“2 


C-j  = 0.5f^  + 0.055 
C2=  0.833 
Pr©=  2.° 

Prt  =0.9 

The  above  equations  are  in 
Favre-averaged  form 
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In  the  fully  turbulent  region  [Y+  = 0(102)],  a compressible 
extension  to  the  high  Reynolds  number  form  of  the  ARSM 
due  to  Galmes  and  Lakshminarayana  (1984)  is  adopted: 


//  //  2.  , 

-P ut  Uj  = - -tijPk  - pkT.j 


where 


+ (P.j  - 2P6(/3)(1-C2) 
* P + pe(Cj-l) 


„ ~ / //  //  //  /A 

“ e*tfPM.  “y  > 


(//  /I 

-pw<  K* 


6bct 


2P 


where  Cj  = 1.5,  C2  = 0.6. 
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In  the  fully  turbulent  region  [Y+  = 0(102)],  a compressible 
extension  to  the  high  Reynolds  number  form  of  the  ARSM 
due  to  Galmes  and  Lakshminarayana  (1984)  is  adopted: 


nu  2 * . j „ 

P uj  = - - pkTi} 


where 


K,P-C2)I2  - (Ptj  - 2P6,/3)(1-C2) 
1 ij 


P + pe(Cj-l) 


Rik  = ^p(eipjpu''u.  - e^jpu-u.) 


//  /A 


/ //  // 
(“P«f 


du . 
; 


//  // 
puy  Z/Jt 


where  Cj  = 1.5,  C2  = 0.6. 
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RKCC  Code 


.Explicit  Multistage,  Conservative,  Compressible 
Formulation 

•Generalized  Coordinates;  2-D,  3-D,  Unsteady; 

— > Turbomachinery  (Rotation  Periodic  B.C.) 

.Coupled  Compressible  Low-Reynolds-Number  K-£ 
Model,  q-co  Model,  ARSM  Model 

•FNS  + Full  Turbulent  KE  Production  Term 

•IRS  for  2-D 

•Finite  Diff  ( Flux  Evaluation  = Cell  Vertex  Finite 
Volume) 

•Characteristic  B.C.’s,  Embedded  H Mesh,  Tip 
Clearance  Topologies  for  Turbomachinery 
Application 

•Eigenvalue  & Local  Velocity  Scaling  of  Artificial 
Dissipation 
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Laser  Measurement  Plane  I II 


Impeller  Inlet  (i  - 15) 


Inlet  Boundary  (i  «=  I ) 


SI 


Diffuser  Exit  (i  * 59) 


Impeller  Exit/ 
Diffuser  Inlet  (i  * 49) 


m 


Figure  a)  Meridional  view  of  59  x 27  x 27  computational  grid  for  Krain  impeller 
computation.  Streamwisc  grid  indexing  and  location  of  experimental  laser  planes  are 
indicated  for  reference,  b)  Cross  stream  views  of  59  x 27  x 27  computational  grid  for 
Krain  impelier  computation. 
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Figure  2 ] ?p  data  acquisition  Plane  IV:  a)  Computed  normalized  relative  helicity 
contours,  b)  Meridional  velocity  profiles.  Experimental  measurements  (symbols), 
computed  values  (solid  line),  c)  Meridional  velocity  contours.  Experimental  contours 
(top),  computed  contours  (bottom). 
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Figure  5 Near  wall  relative  radial  and  crossflow  velocity  profiles  at  90  % chord,  a) 
midspan  on  suction  surface,  b)  midspan  on  pressure  surface,  c)  mid  passage  on  hub. 

Comparison  of  solutions  using  three  turbulence  models  : low  Reynolds  number  k-e  model 

(square),  k-E/ARS  model,  Rjj  = 0 (circle),  k-e/ARS  model,  Rjj  * 0 (triangle).  W3  is  the 
radial  component  (+  towards  tip),  W2  is  the  relative  circumferential  component  (+  towards 
pressure  surface) 
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Experiment  A 


Experiment  B 


Experiment  C 


Figure  6 Schematic  of  experimental  configurations. 
(Note:  z = 0 is  located  .04  chord  inboard  of  the  cascade 
symmetry  plane  [k=nk].  For  the  cases  with  clearance, 
this  corresDonds  to  the  blade  tip  [t/c  = .04]) 


Ff k 


inlet 

plane 


Figure  7 , Views  of  computat 
nomenclature  a)  Detail  of  leading  t 
three-dimensional  grid. 
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Figure  8 . Comparison  of  computed  and  measured  blade 
static  pressure  coefficient  at  three  spanwise  locations. 
Measured  values  (symbols),  Navier-Stokes  (solid  line), 
Euler  (short  dash  line). 
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Figure  12  Comparison  of  computed  and  measured 
spanwise  distribution  of  mass-weighted,  pitchwise 
averaged  outlet  flow  angle.  Measured  values  (symbols), 
Navier-Stokes  (line),  a)  Experiment  B,  t/c  = 0.00. 
b)  Experiment  B.  t/c  = 0.04. 


Figure  13  Comparison  of  computed  and  measured 
spanwise  distribution  of  blade  normal  force  coefficient. 
Measured  values  (symbols),  Navier-Stokes  (line), 
a)  Experiment  B,  t/c  = 0.00.  b)  Experiment  B,  t/c  = 
0.04. 
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Figure  14  Contours  of  flow  angle  at  outlet 
measurement  plane  for  Experiment  C,  t/c  = 0.04. 
a)  Measured,  b)  Predicted,  c)  Comparison  of  computed 
and  measured  spanwise  distribution  of  mass-weighted, 
pitchwisc  averaged  outlet  flow  angle  for  Experiment  C, 
t/c  = 0.04.  Measured  values  (symbols),  Navier-Stokes 


VKI  CASCADE 

Highly  loaded  transonic  turbine  nozzle 
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Figure  15 


Figure  i6  Computational  grid  system  (129*61) 
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Figure  is  The  computed  and  measured  wakes  for 
above  case 

( 0.433  chord  length  downstream) 


Figure  i9  . Mach  no.  contour,  for  above  case 
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Figure  20  Blade  heat  transfer 

(M2is=0.93,Po1=0.915  Bar,  M 1 = 0 

Tu„=1  %,T_ < =403K,  Re, o=  600,000) 
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Figure  21  Computation  Grid  for  a Flat  Plate  Cascade 
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Figure  22  Unsteady  Blade  Surface  Pressure  Jump 
a)  Amplitude,  b)  Phase  Angle 


STATIC  PRESSURE 


77089. 


73490. 


Figure  23  Unsteady  Blade  Pressure  for  a Flat  Plate 
Cascade  ( at  zero  steady  incidence)  with 
an  Incoming  Moving  Wake 
(Runga-Kutta  Explicit  Technique,  M=0.7) 


CONCLUSIONS 

RKCC 


•3D  CODE  WITH  K-e/ARSM  MODEL,  GOOD 
ACCELERATION  PROPERTIES  AND  OPTIMUM 
ARTIFICIAL  DISSIPATION  HAS  BEEN  CODED, 
VALIDATED  AND  USED  EXTENSIVELY  TO  PREDICT 
VISCOUS  FLOW  FIELD  IN  SUPERSONIC  CASCADE, 
SUBSONIC  CASCADES,  LOW  SPEED  COMPRESSOR, 
EVOLUTION  OF  TIP  CLEARANCE  FLOW,  HIGH  SPEED 
CENTRIFUGAL  COMPRESSOR  FLOW  FIELD. 

•3D  CODE  IS  PRESENTLY  EXTENDED  TO  INCLUDE 
TIME  ACCURATE  SOLUTION  ROTOR/STATOR 
INTERACTION  EFFECTS. 

•TURBULENCE  MODEL  SOURCE  TERMS  DO  NOT 
AFFECT  THE  STABILITY  OF  THE  NUMERICAL  SCHEME 
- CONTRARY  TO  THAT  GENERALLY  PERCEIVED. 

•IMPLICIT  TREATMENT  OF  TURBULENCE  SOURCE 
TERM  DID  NOT  IMPROVE  CONVERGENCE  RATE. 

•NO  ADVANTAGE  IS  NUMERICALLY  COUPLING  THE 
GOVERNING  EQUATION  AND  K-e  MODEL. 

•NEAR  WALL  TURBULENCE  PHYSICS  AND  THE 
EFFECT  OF  ROTATION  AND  COMPRESSIBILITY 
PREDICTED  ACCURATELY. 
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•THE  EFFECT  OF  ROTATION  AND  ENDWALL 
SECONDARY  FLOW  IN  THE  NEAR  WAKE,  AND  EFFECT 
OF  ROTATION  ON  3D  BLADE  BOUNDARY  LAYERS  ARE 
CAPTURED  ACCURATELY. 

•THE  SPANWISE  DISTRIBUTION  OF  LOSSES  IS 
PREDICTED  WELL  AWAY  FROM  THE  ENDWALLS,  BUT 
ONLY  FAIR  PREDICTION  IS  ACHIEVED  IN  THE 
ENDWALL  REGION. 

•STAGE  PRESSURE  RISE,  MERIODIONAL  VELOCITIES 
IN  A CENTRIFUGAL  COMPRESSOR  IS  PREDICTED  TO 
ENGINEERING  ACCURACY.  THE  CRITICAL  ELEMENT 
IS  GRID,  FOLLOWED  BY  TURBULENCE  MODELLING. 

•THE  ’ PINCHED"  H GRID  IS  THE  MOST  APPROPRIATE 
GRID  TO  USE  IN  THE  TIP  CLEARANCE  REGION. 

•THE  CODE  ACCURATELY  CAPTURES  SHOCK 
WAVE/BOUNDARY  LAYER  INTERACTION  IN  A 
CASCADE,  WITH  THE  EXCEPTION  OF  REFLECTED 
WAVES  FROM  THE  SUCTION  SURFACE;  SHOCK 
LOSSES  ARE  PREDICTED  WELL. 
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N92~32271 


ABSTRACT 

DEVELOPMENT  OF  A CFD  CODE  FOR  INTERNAL  FLOWS 
IN  LIQUID  FUILED  ENGINES 


SUBMITTED  TO 

WORKSHOP  FOR  COMPUTATIONAL  FLUID  DYNAMIC 
APPLICATIONS  IN  ROCKET  PROPULSION 


To  support  the  design  efforts  of  engines  in  which  liquid  propellants  are  used,  one 
is  often  required  to  analyze  incompressible,  two-dimensional  or  axisymmetric  flows 
within  ducts  and  cavities  with  rotating  walls  and  complicated  geometries.  The  steady- 
state  solution  is  of  interest  in  most  cases. 

This  code  is  intended  to  provide  a tool  for  efficient  CFD  analysis  of  such  flow 
problems,  taking  advantage  of  the  artificial  compressibility  concept  and  the  Beam- 
Warming  numerical  scheme  modified  for  second  order  accurate,  implicit  boundary 
conditions.  These  concepts  ensure  a stable,  robust,  and  accurate  algorithm  due  to 
the  reduced  speed  of  sound  and  the  accuracy  of  the  boundary  conditions. 

The  code  is  dedicated  only  to  two-dimensional  or  axisymmetric  flows  with  or 
without  swirl.  Three-dimensional  computation  is  excluded  to  increase  efficiency  and 
speed. 

This  paper  briefly  presents  the  theory  of  the  code,  as  well  as  several  benchmark 
applications  with  comparison  to  well  known  analytical  solutions.  In  all  these  test 
cases,  the  code  produced  remarkably  accurate  results. 

Youssef  Dakhoul,  Sverdrup  Technology 
620  Discovery  Drive 
Huntsville,  AL  35806 
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TECHNICAL  OBJECTIVES 
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Governing  Equations 


Sverdrup 


Nondimensionalization 
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Pr*t  = P = fluid  density  tn{  = xnt  / vn{  R,  = xnfvn(/  vnt 


Flux  Jacobians 
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Treatment  of  the  Boundary  Conditions 


£ 

tS 

co 

£ 

0 

o 

.2 

<3 

•*3 

M 

CP 

u 

<s 

1-1 

3 

H 

O 

2 

• *«* 

II 

• M 

O 


2 

o 


CP 

-o 

o 

£ 

-£ 

o 

0) 

-*3 

cS 

Ou 

CP 

-4J 

oa 


CO 

b* 

CE 

V 

-£ 

*3 

bO 

£ 


6 

«- 

c2 

_x 

'C 

Is 

6 

_X 

CJ 

_o 

3 

<S 

.2 

TJ 

OP 

BO 

CO 

V 

tH 

a. 

x 

CP 

cp 

£ 

<S 

CJ 

<S 

.B 

-4-3 

CO 

£ 

.2 

'-4—3 

«S 

£ 

a* 

op 


0) 

CO 


CO 

' CP 
CJ 

£ 

-o 

O 


£ 

CP 

> 

CP 


co 

0) 

CJ 


cS 

2 


cp 

-c 


V 

cj 

£ 

■ CP 

-4-J 

60 

‘x 

CP 

cp 

-fl 


(S 

-£ 

-4-3 

-4-3 

o 

Z 


CO 

CP 

bO  co 
-O  cS 
0)  V 


£ 

"<S 

2 


<s 

£ 

O 

bO 

<s 


to 

a 

‘v 

-o 

2 

o 


X 

‘tZ 

ts 

2 

-x 

cj 

o 


cp 

-c 


£ 
o 
0 -A 

o 

^ 3 

*T3  *Bi 
CP  ° 


rj 

"*5 


♦ — * c*  * co 

Crf  o»  c* 

v-o  >-0 


£ 

'o 

o 

CO 

3 


CO 

CP 

.2 

•E 

Is 

2 


CO 

CP 

bO 

■-O 

-4-3 

CO 

V 

* 

-o 

£ 

<S 


3 

a; 

3 


V CJ 

3 * 
H 3 


Q* 

‘O 


o 


o • • • 

• • • 

• • 


CP 

-a 

L- 

0 

1 

"d 

£ 

O 

0 

CP 

CO 

-a 

CP 

-o 

*53 

1 

OP 

£ 

O 

OP 

co 

£ 

<S 

CJ 

CP 

J=> 

u 

<s 

CP 

o. 

a 

<s 

>3 

CP 

-fl 

H 


CO 

CP 

bO 

-O 

CP 

CP 

co 

OP 

_£ 


(S 

-a 

CP 

CO 

£ 

3 

£ 


H 

a 

E 


6* 

MS 


CO 

b* 

cd 

<U 

-c 


co 

Stj  » 

52  2 
| £ S 

£ J3  *3— 
3JD  *3 
v-  £ fl 

© " g 

^ O fa 
0)«jO 
co  TO 
<D 


to 

.2 

‘u 

£ 

CP 

l- 

JP 

te 

the  third  : 
o 8Q\  can 

*6 

b0-4^> 

CP 

cS  CJ 

— 3 

Is 

-.*■£ 

M 

X O 

£ 

CJ 

•g  »- 

J3  a 

3 

«s  w 
2 -c 

MM 

CO 

**■ 

CP 

CO 

k-4 

£ 

O 

U 


cS 

CP 

a 

cS 

co 

CO 


-X  - co 
CJ  > «-4 

-2  I 

•©  z 

CP 

MM  g 

o -O 

Cm  " 

° Is  ^ 

!!  5 2 

S -M3 

MM 

to  aj  u 
£ C £ 

"rf 

£ >3  -fl 

fl  JD  *3 


cS 

i=  VJ» 

72 


▼ 

a o 


OP 


■**  BO  CP 
CP  C 5 
J3  .2 
-5  -fl 

2 £ 
£ cS 
o cr 

*3  CP  J? 
co  *^5 
QJ 

- ° 2 

£ *8  2 


a 

9 

i 

fc 

; 

M 


1318 


last  row  and  for  edges  of  blocked-out  areas  is  required. 
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SLIDE  BLOCK  FLOW 
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SLIDE  BLOCK  FLOU 
PRESSURE  CONTOURS  (Pascal) 
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U Velocity  (m/sec)  at  x/L=0.5 
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Initial  Uelocity  Field  = 0.0 
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Monitor  Point  (51,16 


TRVLOR  UORTICES 

SWIRL  UELOCITY  CONTOURS  tn/sec ) 


■'"’“Oaacaaaaaaaca: 
IQGSaaCCDODDQDCS: 
I !++  + + ■♦■+■♦■  + + + ♦■*-■*■■ 

JUJLlIUJUJUJLiJLJIjJUJLlILlIUJUJLuLJ 

iQQQQDCDDaGsaacQ; 

jooQcacoooDOQccn: 

luasQQcoaoaQQaca; 

IC3aaCM^riCQ3QCvi,TlDOOaCs!~- 

■ to  ro — r-c\cMC'jfMCvinr:ir3; 


, - = = cac  = 
:;:=CCS3 
•*  •+■  + + + + 

u-  L_l  LiJ  UJ  UJ  JJ 

a a a a a a a 
a a a a a e s 
;aaaaaa= 

'^GCMTfflX 
r l~  imn  in  -n 


z 1 GGQQDDDDCDDOQQDDCQ3  ':aaaaaa  = 

o 

u 


1332 


o.6noouE+im 
II.B2000E+0II 
n.nniioF  +nn 


PRESS 


THVLOR  UORTICES 


CO  GO  OB  CO  CO  Cl  Ol  ff)  05  05  05  05  01  01  tC  05 

(2aaaa=aaaacaaaa<-=j: 

->  i i i i i i i i i i i i i i i i 

ujUJUJUUJkJULiJULLllUUJUJUJUU i 

jDDOQODQQDODDODir  = ' 
oaQaeanaanoaaaT  3 
Q-aaaac:naaacaaaa^-~ 
5TnN>-CaDQQCQaQDN3 
Qf-r-i-i O105N.tDLD,TnCSJ<-(T  — 


o:  cncn  oi  m oe  ox  co  x 
zixxcxcxxxx 

i i i i i i i i i i 
_ IUJUUJL-UjJUU 

xxxccexxxx 

zixxxcxxxxa 

xxxaxxxxex 

OXXOCC’-'TSICVT 
' 13  r-  CO  05 


z i i i i i i i i i i i i i i cx  xxxxacxxxx 

a 

u 


1334 


ORIGINAL  PAGE  IS 

OF  POOR  QUALITY 


co  co  03  o 


{enpiso^  gjnssajj  jo  SHH 


1335 


Iteration  (thousands) 


DRIVEN  CAVITY  FLOW 


Initial  Pressure  Field  = 0.0 
Initial  Velocity  Field  = 0.0 


DRIUEN  CflUITV  FLOU 
PARTICLE  TRACES 


337 


I 

o 


in 

I 

O 


CD 

I 

o 


I 

O 


■[ienpis9£[  ajnsssjj  jo  SHH 


CD 


1338 


10- 


TWO  ECCENTRIC  CYLINDERS 
GRID  81  x 31 


14 

0) 

•o 

c 


S 

o 

M 


1339 


Initial  Velocity  Field 
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HTD  PREBURNER  DISCHARGE  CHUITV 
Relative  Pressure  Contours  (psi) 
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RTD  PREBURNER  DISCHARGE  CRUITV 
Swirl  Uelocity  Contours  (fps) 
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Development  of  the  KIVA-II  CFD  Code  for  Rocket  Propulsion  Applications 

by  Robert  V.  Shannon  Jr. 
and 

Alvin  L.  Murray 

Aerotherm  Corporation 
1500  Perimeter  Parkway,  Suite  225 
Huntsville,  AL  35806 

Due  to  the  existence  of  liquid,  solid,  and  hybrid  rocket  motors  a need  exists  for  a fluid 
dynamics  code  which  can  solve  for  both  the  Navier-Stokes  equations  in  multi-dimensions 
as  well  for  liquid  and  solid  particle  motion  within  the  gas  flow  domain.  This  type  of  CFD 
code  must  couple  the  gas  and  particle  motion  so  that  the  effects  of  one  upon  the  other 
can  begin  to  be  understood.  Disciplines  such  as  the  evaluation  of  solid  motor  performance 
as  well  as  nozzle  erosion  predictions  for  solid  motors  have  a great  need  for  the  type  of 
information  that  this  type  of  computer  code  can  provide.  In  addition,  it  is  difficult  to 
accurately  simulate  liquid  motor  combustion  chamber  flows  without  solving  for  the  liquid 
oxidizer  droplet  motion,  breakup,  and  evaporation  coupled  with  the  reacting  gas  flow. 

The  KIVA-II  code,  originally  developed  at  Los  Alamos  National  Laboratories  to  solve  fluid 
dynamics  problems  in  internal  combustion  engines,  has  been  developed  to  solve  rocket 
propulsion  type  flows.  This  work  was  supported  by  the  NASA  Solid  Propulsion  Integrity 
Program.  The  objective  of  the  work  performed  was  to  develop  this  CFD  code  so  that  both 
liquid  and  solid  particle  motion  could  be  simulated  for  arbitrary  geometry  and  for  high 
speed  as  well  as  low  speed  reacting  flows. 

Modification  of  the  KIVA-II  gas  flow  algorithm  involved:  incorporating  independently  speci- 
fiable supersonic  and  subsonic  inflows  and  outflows,  symmetric  as  well  as  periodic  boundary 
conditions,  the  capability  to  use  generalized  single  or  multi-specie  thermodynamic  data  and 
transport  coefficients  and  allowing  the  user  to  specify  arbitrary  wall  temperature/  heat 
flux  distributions.  Major  modifications  to  the  algorithms  involving  both  convection  and 
diffusion  were  successfully  performed  and  results  have  been  compared  with  some  available 
experimental  data  on  nozzle  flows  to  verify  these  modifications. 

Modification  of  the  algorithms  governing  particle  motion  involved:  incorporation  of  multi- 
specie particle  types,  generalized  liquid  property  thermodynamic  data,  options  which  pre- 
vent or  allow  particle  evaporation,  a generalized  particle  injection  algorithm,  simulation  of 
particle  elastic  or  inelastic  collisions  with  solid  boundaries,  periodic  or  symmetric  particle 
boundary  conditions  at  boundaries  which  are  not  inflow,  outflow,  or  wall  boundaries,  and 
incorporation  of  a particle  drag  model  which  allows  for  both  Reynolds  and  Mach  number 
effects  on  particle  drag  coefficients. 

This  new  CFD  code  has  been  shown  to  successfully  solve  rocket  propulsion  flows  as  well 
as  rocket  propulsion  flows  with  entrained  particles  for  several  different  rocket  nozzles,  sub- 
merged and  otherwise.  Verification  of  this  new  CFD  code  continues  by  comparing  results 
to  experimental  data  as  well  as  to  predictions  of  other  CFD  codes.  This  program  is  proving 
to  be  a valuable  tool  in  the  prediction  of  rocket  propulsion  flows. 
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Development  of  the  KIVA-II  CFD  Code 
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Aerotherm  Corporation 

A Subsidiary  of  OynCorp 


Development  of  the  KIVA-II  CFD  Code 
for  Rocket  Propulsion  Applications 
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Aerotherm  Corporation 

A Subsidiary  of  DynCorp 


Original  KIVA-II  Code  - General  Information 
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VELOCITY  VECTORS  FOR  PLANAR  NOZZLE 
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TWO  INCH  MOTOR  NOZZLE  TEST  CASE 
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MNASA  MOTOR  NOZZLE  VELOCITY  VECTORS 


MNASA  MOTOR  TEMPERATURE  FIELD 
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PARCEL  VELOCITY  CHRONOLOGY 
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PARCEL  VELOCITY  CHRONOLOGY 


PARCEL  VELOCITY  CHRONOLOGY 


PARCEL  VELOCITY  CHRONOLOGY 


A new  finite  volume  2D/3D  Navier-Stokes  CFD  code  has  been 
developed  by  using  the  KIVA-II  code  as  the  starting  point. 
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N 92 “32273 

A Computational  Design  System  for  Rapid  CFD  Analysis 

E.P.  Ascoli,  S.L.  Barson,  M.E.  DeCroix,  and  M.M.  Sindir 
CFD  Technology  Center 
Rockwell  International,  Rocketdyne  Division 
Canoga  Park,  California 

Effective  application  of  computational  fluid  dynamics  (CFD)  in  the  engineering 
environment  requires  that  key  design  and  analysis  tools  be  integrated  to  the  greatest 
extent  possible.  A computational  design  system  (CDS)  is  described  in  which  these  tools 
are  integrated  in  a modular  fashion.  This  CDS  ties  together  four  key  areas  of 
computational  analysis;  description  of  geometry,  grid  generation,  computational  codes, 
and  postprocessing.  Common  input  and  output  formats  and  necessary  translators  are 
established  to  facilitate  data  transfer  between  the  four  key  areas  and  to  enhance  system 
modularity.  Advances  made  in  three  key  areas  are  described. 

Significant  progress  has  been  made  toward  integration  of  geometry  definition  systems 
with  CFD  grid  generation  tools.  Geometric  data  have  successfully  been  passed  from  the 
Catia  CAD  and  Patran  CAE  systems  to  the  Rockwell  Automated  Grid  Generation  System 
(RAGGS).  The  IGES  standard  was  employed  in  each  case.  The  CFD  Pump  Consortium 
impeller  geometry  is  used  to  illustrate  that  a reasonable  level  of  integration  is  achievable 
for  complex  geometries. 

While  many  CFD  grid  generators  are  available  in  the  public  domain,  their  capabilities  vary 
widely.  An  effort  is  underway  to  systematically  review  available  codes,  identify  those  most 
applicable,  and  integrate  them  into  the  CDS  described.  Seven  grid  generation  codes  are 
currently  being  reviewed  according  to  previously  defined  evaluation  criteria  described 
herein. 

Postprocessing  tools  are  being  developed,  reviewed,  and  integrated  into  the  CDS  as 
well.  Engineering  data  extraction  tools  are  being  developed  to  be  completely  consistent 
with  existing  code  methodologies  and  to  remain  completely  modular.  Tools  for  data 
visualization  have  advanced  noticeably  over  the  last  few  years.  These  are  being 
reviewed  and  integrated  into  the  CDS. 

Integration  of  improved  CFD  analysis  tools  through  integration  with  the  Rocketdyne  CDS 
has  made  a significant  positive  impact  in  the  use  of  CFD  for  engineering  design  problems. 
Complex  geometries  are  now  analyzed  on  a frequent  basis  and  with  far  greater  ease. 
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INTEGRATION  OF  CFD  TOOLS  NEEDED  FOR  EFFICIENT 
APPLICATION  IN  ENGINEERING  ENVIRONMENT 
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ROCKETDYNE  ADVANCED  COMPUTATIONAL 
ENGINEERING  SYSTEM  (RACES) 
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ENG.  DATA  DATA  SOLUTION 

EXTRACTION  VISUAL.  DATABASE 

Rockwell  International  

Rocketdyne  Division  — - Translators  • - Common  Translator  I/O  Format 


CAD  / CAE  INTERFACES  ESTABLISHED 
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COMPUTATIONAL  DESIGN  SYSTEM  SIGNIFICANTLY 
ENHANCES  CFD  ANALYSIS  CAPABILITY 
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OPTIMUM  DESIGN  OF  NINETY  DEGREE  BENDS 


Vijay  Modi 

Department  of  Mechanical  Engineering 
Columbia  University 


Abstract 


An  algorithm  for  the  optimum  design  of  an  internal  flow  component  to 
obtain  the  maximum  pressure  rise  is  presented.  Maximum  pressure  rise  in  a 
duct  with  simultaneous  turning  and  diffusion  is  shown  to  be  related  the 
control  of  flow  separation  on  the  passage  walls.  Such  a flow  is  usually 
associated  with  downstream  conditions  that  are  desirable  in  turbomachinery 
and  propulsion  applications  to  ensure  low  loss  and  stable  performance.  The 
algorithm  requires  the  solution  of  an  "adjoint"  problem  in  addition  to  the 
"direct"  equations  governing  the  flow  in  a body,  which  in  the  present 
analysis  are  assumed  to  be  the  laminar  Navier-Stokes  equations.  Earlier 
studies  have  usually  addressed  such  problems  for  the  case  of  inviscid 
and/or  irrotational  flow.  These  assumptions  may  not  be  valid  in  flows  that 
undergo  sharp  turning  resulting  in  strong  secondary  flows  and  possibly 
separating  and  recirculating  regions.  The  theoretical  framework  and 
computational  algorithms  presented  in  this  study  are  for  the  steady  Navier- 
Stokes  equations. 

A novel  procedure  is  developed  for  the  numerical  solution  of  the 
adjoint  equations.  This  procedure  is  coupled  with  a direct  solver  in  a 
design  iteration  loop,  that  provides  a new  shape  with  a higher  pressure 
rise.  This  procedure  is  first  validated  for  the  design  of  optimum  plane 
diffusers  in  two-dimensional  flow.  The  direct  Navier-Stokes  and  the 
"adjoint"  equations  are  solved  using  a finite  volume  formulation  for 
spatial  discretization  in  an  artificial  compressibility  framework.  The 
discretized  equations  are  integrated  using  explicit  Runge-Kutta  time  steps 
to  obtain  steady-state  solutions.  It  is  found  that  the  computational  work 
required  to  solve  the  "adjoint"  problem  is  of  the  same  order  as  that 
required  to  solve  the  direct  problem.  It  is  also  found  that  the  procedure 
converges  within  about  ten  iterations,  and  in  addition,  the  number  of 
design  iterations  are  not  sensitive  to  the  grid  used  for  the  calculations. 
This  is  a significant  computational  advantage  over  heuristic  design 
procedures  based  on  point  by  point  sensitivity  analysis  where  the  work 
increases  with  the  refinement  of  the  grid. 

A simplified  version  of  the  above  approach  is  then  utilized  to  design 
ninety  degree  diffusing  bends.  The  bend  inlet  is  square  with  intermediate 
and  exit  cross-sections  constrained  to  be  rectangular.  The  location  of  bend 
walls  is  then  determined  in  order  to  obtain  the  maximum  pressure  rise 
through  the  bend.  Calculations  were  carried  out  for  a mean  radius  ratio  at 
inlet  of  2.5  and  Reynolds  numbers  varying  from  10f)  to  500.  While  at  this 
stage  laminar  flow  is  assumed  it  is  shown  that  a similar  approach  can  be 
conceived  for  turbulent  flow’s. 
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MOTIVATION 


# How  to  shape  internal  flow  passages 

# Combined  turning  and  diffusing  flow 

# Maximize  pressure  rise 

# Can  not  assume  inviscid  or  2-D  flow 

OBJECTIVES 

# Develop  theoretical  framework  — laminar  3-D 

# Develop  Navier-Stokes  and  Adjoint  solvers 

# Validate  Navier-Stokes  solver 

(Laminar  90  degree  bend,  Taylor  et  al.  1982) 

# Validate  Optimization  Approach  on 
2-D  straight  diffusers 

# Apply  to  the  design  of  ninety  degree  bends 
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'U'i,!  — 0 

UjUiJ  = -p*i  + UUiJj  , p*  =p/p 

# No  slip  BC  on  Tm 

# Dirichlet  BC  for  at  Tj  and  To 


p*UiUids  + 

1400 


p*Uin{ds 


J (r*Me ) ~ J(Xm)  — + 0(e2) 


# Our  “Adjoint”  equation: 


Vzi, jj  Zjti)  r,i  0 


in  fi 


m = - & - ui ) 


q = | (r  - P*  + (l/2)ui  - 2«ywy) 


# Pironneau’s  “Adjoint”  equation: 

Wii  = 0 


vwi}jj  + UjWi'j  + WjUj'i  - qti  = —UjUij 


U);  = 0 


in  Q 
in  fi 
on  r 


8J 


= V [ p{s)  ( 


du{\  ( dui  dwi 


1 l 1-  2 - 1 ds 

vM  \dnj\dn  dn 


dw 

dt 


= w (const)  + • • • 
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NAVIER-STOKES  SOLVER 


pt  = -P2uiti 
'U'ij  “h  UjUij  — 

a)  Artificial  Compressibility 

b)  Runge-Kutta  Time  Integration 

c)  Finite  Volume  Discretization 

d)  Artificial  Dissipation 

e)  Local  Time  Stepping 

e)  Implicit  Residual  Smoothing 
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Figure  4:  Geometry  of  a circular  bend  with  square  cross  section. 


Figure  5:  Streamwise  velocities  in  a bend:  a)  0 = 30  degrees,  b)  © — 60  degrees,  c)  0 77.5  degrees. 
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ADJOINT  EQUATION  SOLVER 


MESH  GENERATION 
-#■  Thompson  et  al. 

&££  %TJT]  0 

vet  + yVv  = 0 
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Boundary  Conditions 

# No-slip  be  on  the  walls 

# Zero  normal  derivatives  at  exit 

# Streamwise  velocity  component  is  specified  at 
entrance 

# Zero  normal  derivatives  for  remaining  velocities 
at  entrance 

# Typical  bc’s  at  symettry  planes 

=#=  Zero  second  derivatives  for  pressure  (a  compu- 
tational be) 


Diffuser  Profile  History 


Wall  Shear  Stress 


Skin  Friction  History 


Re=200,  61  by  31  grid 

o : N=1 

n : N=2 
A : N=5 
* : N=10 
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Pressure  Rise 


Pressure  Rise  History 


Re=200,  61  by  31  grid 

O : Area-averaged  pressure  rise 
A : Flow-averaged  pressure  rise 
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Pres 


Actual  Pressure  Rise 


Ideal  Pressure  Rise 


Reynolds  Number 


O • Optimal  diffusers 
A : Straight  diverging  diffusers 


2 


Grid  Study 


T 

2 


Re=200 

O : 31  by  16 
n : 61  by  31 
* : 121  by  61 


The  error  in  the  location  of  the  optimal  diffuser 
profile  corresponding  to  a 2 percent  error  in  the 
total  pressure  rise.  The  optimum  shape  lies  between 
the  high  and  low  y-values  shown  in  the  graph 
(Re=500,  grid  size  is  61  by  21,  and  L/Wi  = 3). 
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Figure  3: 


Sketch  of  the  inlet  header 


A representative  section  of  a three-dimensional  diffuser.  Flow  enters  at  upstream  boundary  Tj 
and  exits  at  downstream  boundary  To.  The  walls  to  be  shaped  are  TM. 
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ISSUES 


# Geometry  Constraints 

In  general:  Move  walls  by  ep(s)  along  the 

normal  direction,  everywhere 

New  shape  may  not  satisfy 

— specified  mean  passage  location 

— specified  cross-sectional  shape 

— overall  system  geometry 

# Present  Work 

- No  correction  on  side  walls  (z  = 0,  z = 2rmax) 

- Apply  mid-plane  (z  = zm ax/2)  correction 
to  all  £ locations 

- Hence  all  cross-sections  are  rectangular 

# Laminar  flow  results  (Re  < 500) 


1417 


Governing  Equations: 


ujuijj  = ”"P,t  "b  l/ui,jj 


Design  Objective: 

Maximize  Static  Pressure  Rise 


Cabuk  and  Modi,  1990 


= 6 


1418 


• • 


DESIGN  ALGORITHM 


1:  Choose  an  initial  shape. 

Generate  the  computational  grid. 

3:  Solve  the  N-S  equations. 

4:  Compute  shear  stress  on  the  walls. 

5:  Compare  wall  shear  stress  to  target  dis- 
tribution and  determine  the  amount  of 
boundary  movement  p(s). 

6:  Update  the  shape. 

7:  Go  to  step  (2) 


1419 


Normalized  wall 
shear  stress 


Wall  shear  stress  along  the  outer  wall, 
Re=100 

O • Optimum  diffusing  bend 
A : Elliptic  diffusing  bend 
Dashed  Curve  : Target  distribution 
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Normalized  wall 
shear  stress 


1.4 


Arclength  along  the  bend 


Wall  shear  stress  along  the  inner  wall, 
Re=100 

O • Optimum  diffusing  bend 
A : Elliptic  diffusing  bend 
Dashed  Curve  : Target  distribution 
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Optimum  Shapes 
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Performance  = f(shape) 

^ Possible  Applications: 

- 90  Degree  Bend 

- Turn  Around  Ducts 
— Transition  Ducts 

- S-shaped  Ducts 

— Straight  or  Curved  Diffusers 
— Turbine  Blades 

- Engine  Inlets 

- Turning  Vanes 
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CONCLUSIONS 


Theory 

# Theoretical  Framework  for  Design 
with  Navier-Stokes  equations 

# Determine  p(s)  from  Direct-f- Adjoint 
or  from  Direct  alone 

Computational 

# Direct  and  Adjoint  Solvers  Validated 
for  Plane  Diffusers 

# Design  of  90  degree  Bend  with  Specified 
Cross-section,  Max.  Ap 

# Number  of  Design  Cycles  < 10 

# “Flow”  Interpretation  of  Adjoint  Problem 

Future  Plan 

# Apply  to  3-D  turbulent  flow 

# Specify  mean  line,  vary  cross-section 

# Other  objectives  : Min.  Distortion 
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A Multi-Domain  Method  for  Subsonic  Viscous  Flows 
Daniel  C.  Chan  and  Munir  M.  Sindir 
CFD  Technology  Center 

Rocketdyne  Division,  Rockwell  International  Corporation 

Canoga  Park,  California 

We  have  developed  a Schwarz  type  domain  decomposition 
method  for  a pressure  base,  two-  and  three-dimensional  Navier- 
Stokes  solver.  This  technique  allows  one  to  partition  a flow  path, 
which  can  be  characterized  by  complex  geometry  and/or  complicated 
flow  physics,  into  smaller  sub-domains  according  to  the  local 
geometric  simplicity  or  estimated  flow  scales.  We  can,  then,  sweep 
the  sub-domains  in  some  order  and  solve  the  Navier-Stokes 
equations  using  as  boundary  conditions,  along  the  domain  interfaces, 
the  Dirichlet  conditions  which  are  taken  from  the  most  recent  update 
of  the  solution  in  the  adjacent  neighboring  domains.  With  this 
technique,  one  can  minimize  the  adverse  effects  caused  by  grid 
skewness  and  the  stiffness  problem  caused  by  disparate  flow  scales. 

This  code  has  been  successfully  applied  to  complicated 
engineering  problems  and  the  results  are  presented  as  separate 
papers  in  this  conference.  Here,  we  report  the  results  of  a few 
fundamental  flow  cases  to  demonstrate  that  a judicious  use  of  the 
multi-domain  method  can  offer  a significant  convergence  acceleration 
over  the  traditional  one-domain  method.  This  method  can  be 
extended  to  exploit  the  architecture  of  a parallel  computer  to  further 
improve  the  speed. 
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A MULTI-DOMAIN  METHOD  FOR  SUBSONIC 

VISCOUS  FLOWS 
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MULTI-DOMAIN  PATCHING  ALGORITHM 
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FLOW  OVER  A HALF  CIRCLE 

SINGLE  AND  FOUR  DOMAIN  GRID  TOPOLOGY 
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FLOW  OVER  A HALF  CIRCLE 

SINGLE  AND  FOUR  DOMAIN  FLOW  SOLUTION 
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GREEN'S  FUNCTION  MULTI-DOMAIN  ALGORITHM 

ONE-DIMENSIONAL  HELMHOLTZ  EQUATION  FORMULATION 

PULICANI,  1988 
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GREEN'S  FUNCTION  MULTI-DOMAIN  ALGORITHM 
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ONE-DIMENSIONAL  HELMHOLTZ  EQUATION 
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ONE-DIMENSIONAL  HELMHOLTZ  EQUATION 
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2-D  LAMINAR  BACKWARD  FACING  STEP 

PREDICTED  STREAMLINE  AT  DIFFERENT  REYNOLDS  NUMBERS 
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2-D  LAMINAR  BACKWARD  FACING  STEP 

REATTACHMENT  LENGTHS  AT  DIFFERENT  REYNOLDS  NUMBERS 
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2-D  CAVITY  WITH  AN  ORIFICE 
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MODELING  OF  COMPLEX  FLOW  PATHS 

• GEOMETRIC  COMPLEXITY 
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LARGE  EDDY  SIMULATION  OF  COMPRESSIBLE 
TURBULENT  CHANNEL  FLOWS* 

Robert  A.  Beddini  and  Jeffrey  P.  Ridder 
Department  of  Aeronautical  and  Astronautical  Engineering,  MC  236 
University  of  Illinois  at  Urbana-Champaign 
Urbana,  Illinois  61801 


Statement  of  Problem:  The  development  of  turbulence  within  rocket  propulsion  chamber  flows 
remains  a difficult  problem  to  predict.  Within  solid  propellant  rockets,  for  example,  the  flow  can 
exhibit  multiple  regions  of  transition  to  turbulence,  and  is  susceptible  to  various  modes  of 
aeroacoustic  interaction,  potentially  associated  with  instability  of  the  combustion/flowfield 
process. 

Objective:  To  formulate,  develop  and  validate  a large  eddy  simulation  (LES)  method  for 
compressible  channel  flows.  Additionally,  to  assess  the  potential  and  limitations  of  the  method 
with  regard  to  predicting  flows  of  interest  in  realistic  systems. 

Approach:  The  LES  method  separates  the  resolvable  scale  motions  from  the  unresolvable  (subgrid) 
scales  by  applying  a spatial  filter  to  the  compressible  Navier-Stokes  equations.  The  subgrid-scale 
Reynolds  terms  are  modeled  using  a compressible  extension  of  an  existing  incompressible  model 
for  wall  bounded  flows.  The  equations  are  solved  numerically  using  a modified  four-step  Runge- 
Kutta  procedure  in  time  and  second  or  fourth-order  finite  differences  in  space. 

Results  and  Conclusions^:  The  method  has  been  validated  by  simulating  a low  Reynolds  number 
(Reb  « 5400),  low  Mach  number  (Mc  » 0.3)  turbulent  Poiseuille  flow.  Various  statistical 
comparisons  are  made  with  incompressible  experimental  and  direct  simulation  data  at  similar 
Reynolds  numbers,  including  higher-order  statistics  and  spatial  correlations.  The  results 
compare  favorably  with  the  incompressible  data. 

A high  subsonic  Mach  number  (Mc  = 0.3)  turbulent  Poiseuille  flow  is  also  simulated  for 
comparison  with  the  low  Mach  number  results  at  nominally  constant  Reynolds  number.  The 
mean  velocity  profile  is  seen  to  depart  from  the  low  Mach  number  profile,  corresponding  to  an 
expected  dependence  of  the  mean  density  and  temperature  profiles  on  Mach  number.  The 
turbulence  velocity  statistics  are  found  to  be  reasonably  independent  of  Mach  number.  Pressure 
fluctuation  statistics  are  also  found  to  scale  with  the  wall  shear  stress  independently  of  Mach 
number,  although  normalized  density  and  temperature  fluctuations  increase  substantially  with 
Mach  number.  The  density  and  temperature  fluctuations,  although  small  in  magnitude,  are 
observed  not  to  be  isobarically  related. 

The  current  simulations  have  validated  the  algorithm  in  the  incompressible  limit  and 
have  demonstrated  the  ability  of  the  method  to  simulate  high  subsonic  Mach  number  flows.  The 
principal  impediment  in  application  of  the  method  to  practical,  high  Reynolds  number  chamber 
flow  problems  is  the  large  CPU  time  requirement  of  the  calculations.  At  present,  this  bottleneck  is 
caused  in  large  part  by  resolution  requirements  for  the  turbulence-producing  "streaks"  in  the 
viscous  wall  layer.  Improvements  in  the  subgrid-scale  modeling  could  greatly  reduce  CPU 
requirements,  leading  to  more  rapid  engineering  use. 


* Research  supported  by  NASA  Marshall  Space  Flight  Center  under  grant  NGT-  50363. 

* Ridder,  J.P.:  Large  Eddy  Simulation  Of  Compressible  Channel  Flow,  Ph.D.  Thesis,  U.  of 
Illinois  at  Urbana-Champaign,  January,  1992. 
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APPLICATION  CONDITIONS 
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Turbulent  Channel  Flow 
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RESULTS  (HIGH  Mach  No.) 
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MEAN  VELOCITY  PROFILE  COMPARISON 
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RESOLVED  REYNOLDS  SHEAR  STRESS 

COMPARISON 
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RESOLVED  RMS  PRESSURE  AND 
TEMPERATURE  FLUCTUATIONS 
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N92-32277 

Treating  Convection  in  Sequential  Solvers 


Wei  Shyy 
Siddharth  Thakur 


Department  of  Aerospace  Engineering, 
Mechanics  and  Engineering  Science 
University  of  Florida,  Gainesville,  FL 


The  treatment  of  the  convection  terms  in  the  sequential 
solver,  a standard  procedure  found  in  virtually  all 
pressure  based  algorithms,  to  compute  the  flow  problems 
with  sharp  gradients  and  source  terms  is  investigated. 
Both  scalar  model  problem  and  one-dimensional  gas 
dynamics  equations  have  been  used  to  study  the  various 
issues  involved.  Different  approaches  including  the  use 
of  nonlinear  filtering  technique  and  the  adoption  of  TVD 
type  schemes  have  been  investigated.  Special  treatments 
of  the  source  terms  such  as  pressure  gradients  and  heat 
release  have  also  been  devised,  yielding  insight  and 
improved  accuracy  of  the  numerical  procedure  adopted. 
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A Proposed  Hierarchy  of  Test  Problems 


1.  Nonlinear  scalar  wave  equation  to  study  capturing  of  sharp  gradient 


2.  1-D  gas  dynamics  in  a tube  

(a)  to  study  coupling  among  u,  v,  p and  p 


(b)  to  study  interaction  with  B.C.s  and  among  waves 


open 

end 


closed 

end 


p = 3.5277 

p = 0.571 

p = 0.445 

p =0.5 

u = 0 

closed 

end 


Diaphragm 


(c)  to  study  formation  and  propagation  of  acoustic  and  entropy  waves 


closed 

end 
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3.  1-D  Combusting  Flow  in  a Duct 

to  study  effect  of  chemical  heat  release  and  its  impact 

on  acoustic  and  entropy  waves. 

to  study  propagation  and  interaction  of  nonlinear 

waves  in  pressure,  thermal  and  convective  fields. 

to  understand  the  longitudinal  combustion  instability 

and  to  devise  active  control  strategy 


^Gutter 

Location 


4.  Multidimensional  Problems 
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| For  a numerical  scheme  two  features  control  its  performance  : 

• Amplification  factor ^ numerical  damping 

• Phase  angle  numerical  dispersion 


■ Problems  of  convection  treatment : 

• First  order  upwind  scheme ■-  excessive  damping 

dispersion  problem  suppressed 

• Higher  order  schemes  — ► no  excessive  damping 

dispersion  problem  appears 


■ Two  approaches  investigated 

• Nonlinear  filtering 

• TVD  type  approach  with  source  term  treatment  and 
artificial  compression 
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Backward  Euler  time  stepping  scheme  : 


O')  C = 0.25 
(I.)  C = 0.5 


•(to— ' 


-i  ■ Exact 
solution 


3 i t 1 » 10 


2ic/f)(  = \vavclciiv||i/A.\) 


o \ : j 


2ir/fl(  = iv;tvclciigtli/Ax) 


first-order  upwind  scheme 


•>; (n) 


(to 


(n)  C = 0.2S 
(10  C = 0.5 


2 it/ P(  = wnvcl«:.iBlli/A*) 

2it/|l(  = wnvcIcng(Ii/Ax) 

second-order  central  difference  scheme 



0 

\ 

-O.  3 

{n) 

,(n) 

0 

(t>) 

1 *l 

0>) 

/ 

(:i)  C = 0.25  ' 

-1.3 

(«) 

00  C = 0.5 

('») 

) 1 1 ( 1 1 » 1 
2it/fl(  = wavelength/ Ax) 

( 

0 

» 1 l > * 3 « t 

2re/[)(  = wavelength/ Ax) 

second-order  upwind  scheme 

1 


Point  : A highly  dispersive  scheme  may  become  a good  one  if  dispersive 
problems  in  high  wave  number  can  be  fixed;  or  better  yet,  to  make 
constructive  use  of  these  wiggles 

A possibility  : Nonlinear  Filtering 
Key  elements  : 

* maintains  conservation  laws 

* utilizes  standard  schemes  as  basis 

* attempts  to  eliminate  wiggles  a posteriori  ( a geometric  approach) 

* effective  only  for  short  wavelength  oscillations  (2A  and  4A)  & hence  can 
check  the  filtering  effectiveness  via  grid  refinement 

Define 

Energy  Content 


Area  Content 
— Goal  : 

* minimize  E 


* maintain  A 
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0.6  - 
0.4  - 
0.2  - 

120  160  200  0 


Effect  of  filtering  on  a 1-D  compressible 
Re  = 104 


40  SO 

x 

flow  solution, 


time  step  = 600 


m 

F = mVP  + P 
(E  + p)m/  p 


Can  also  write 


Speed  of  sound 


The  eigenvalues  of  Jacobian  matrix  are 

{a1, a2, a3)  = ( u-c , u,  u + c ) 
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Sequential  Approach  with  Coordinated  Chracteristics 


ap_  + a[p  u]  = 0 

dt  dx 


dm 

d [mu] 

_dp 

dt 

dx 

dx 

dE 

+ d[Eu]  _ _ 

d(pu) 

dt 

dx 

dx 

Pressure  terms  : source  terms 

The  local  characteristic  speed  : convection  speed 

ahV  2 = j(ui  + Um) 


(same  for  all  equations) 
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Special  Source  Term  Treatment 

Conservation  law  with  a source  term  Y(u) 

ut  + /M,  = 'I'M 

Examples: 

Method  I - MacCormack’s  explicit  predictor-corrector  method 

Method  II  : Operator  splitting  (Strang’s  time-splitting) 

t/"’1  = Sv(k/2)  Sj(k)  St(jfc/2)  V" 

where  Sf  represents  the  numerical  solution  operator  for  the  system 

uc  + /(O*  = 0 

and  is  the  numerical  solution  operator  for  the  ODE 

u(  = tO) 
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a)  First-order  upwind  and  second-order  Lax-Wendroff  schemes. 


(b)  Harten’s  TVD  scheme  (6  = 0). 


Density  profiles  using  the  simultaneous  solution  approach  using  different  schemes 


boundary  boundary 


Conclusions 


For  sequential  solvers,  coordination  of  propagation  speed 
among  equations  requires  extra  care. 


Can  apply  modern  TVD  type  schemes  with  source  term 
treatment  to  improve  accuracy. 


Can  utilize  nonlinear  filtering  techniques  to  eliminate 
dispersion  problems. 


Several  test  problems  have  been  investigated; 
results  show  promise. 
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